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ABSTRACT 

We analyze the physical conditions in the interstellar gas of 11 actively star-forming galaxies at z^2, 
based on integral- field spectroscopy from the ESO-VLT and HST/NICMOS imaging. We concentrate 
on the high Ha surface brightnesses, large line widths, line ratios and the clumpy nature of these 
galaxies. We show that photoionization calculations and emission line diagnostics imply gas pressures 
and densities that are similar to the most intense nearby star-forming regions at z=0 but over much 
larger scales (10-20 kpc). A relationship between surface brightness and velocity dispersion can be 
explained through simple energy injection arguments and a scaling set by nearby galaxies with no free 
parameters. The high velocity dispersions are a natural consequence of intense star formation thus 
regions of high velocity dispersion are not evidence for mass concentrations such as bulges or rings. 
External mechanisms like cosmological gas accretion generally do not have enough energy to sustain 
the high velocity dispersions. In some cases, the high pressures and low gas metallicites may make 
it difficult to robustly distinguish between AGN ionization cones and star formation, as we show for 
BzK-15504 at z=2.38. We construct a picture where the early stages of galaxy evolution are driven by 
self-gravity which powers strong turbulence until the velocity dispersion is high. Then massive, dense, 
gas-rich clumps collapse, triggering star formation with high efficiencies and intensities as observed. 
At this stage, the intense star formation is likely self-regulated by the mechanical energy output of 
massive stars. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: kinematics and dynamics 
— infrared: galaxies 



1. INTRODUCTION 

Elucidating the physical processes that regulate global 
star formation is one of the keys to understanding galaxy 
evolution. The power-law relation, over several orders of 
magnitude, between the warm HI and cold molecular gas 
surface density and star formation intensity ("Schmidt- 
Kennicutt law") is telling us that there must be under- 
lying physical processes that control and regulate star 
formation (e.g.. iKennicutt et al.l 2007). However, with a 
myriad of possible mechanisms for regulating star forma- 
tion over large scales - cloud formation and destruction, 
mechanical energy output from stars and AGN, spiral 
density waves, turbulence induced by gravity and me- 
chanical energy, magnetic fields, mixing layers, and many 
others - it is challenging to distill a unifying explana- 
tion over many orders of magnitudes in gas density. One 
possible way of advancing our understanding of the pro- 
cesses that regulate global star formation is by studying 
the most extreme star-for ming galaxies locally and at 
cosmological distances (e.g.. lLehnert k, Heckman f l996a: 

1 Data obtained as part of Programme IDs 075.A-0466, 076. A- 
0527, 077.A-0576, 078.A-0600, and 079.A-0341 at the ESO-VLT 
and based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele- 
scope Science Institute, which is operated by the association of 
universities for research in astronomy, inc., under NASA contract 
NAS 5-26555. 

2 Present Address: Institut dAstrophysique Spatiale, UMR 
8617, CNRS, Universite Paris-Sud, Batiment 121, F-91405 Orsay 
Cedex, France 



iShaplev et al~l 12001 [Verma et al J 12001 . 

Thanks to recent technological improvements, we are 
now able to study the basic physical processes driving 
galaxy assembly directly and over large ranges of cosmic 
time. In particular, the number of galaxies at redshifts 
of 2— 3 with detailed integral-field spectroscopic studies 
in the rest-frame optical is still small, less than a few 
dozen, but growing rapidly. These observations allow us 
to trace spatially resolved emission- line properties, and 
to investigate galaxy kinematics and the physical condi- 
tions of the ionized gas. Up to now, most studies have 
focused on the kinematic properties of the galaxies and 
less so on the properties of the emission lines themselves 
beyond their relative velocities and widths. Recombina- 
tion line fluxes were used to estimate star-formation rates 
using si mple prescriptions d e veloped for low-redshift 
galax i es (j Gcnzcl ct al. I [20061: IForster Schreiber et al. I 
20061: Ivan Starkenburg et al. I [2(jM lLaw et all 120071 : 



Wrig ht et al. Il2007t: INesvadba et al. Il2008al) . 

With these estimates and assumptions, far-reaching 
conclusions have been made regarding the modes of 
galaxy assembly and the drivers of star formation 
in the early Universe. Many of these studies (e.g., 
iGenzel et al~ll2006t IForster Schreiber et al. 1120061 ) favor 
a scenario where high-redshift galaxies exhibit gaseous 
disks of ^10-20 kpc in size, which, unlike galaxies at 
low redshift, are dynamically hot with large gas veloc- 
ity dispers ions, a, relative to th eir bulk velocities, v 
(see also INesvadba et al.l l2006aT) . The ratios of ran- 



2 



Lehnert et al. 



dom to large scale velocity shear at high redshift are 
of order v/tr ~few, compared to v/ cr > 10 at low red- 
shift. To explain the s e obse r vations. iKeres et al. I (2005, 
f2008h : lOcvirk et all (|200ll : iDekel et al. I (|2009bh pro- 
posed a scenario where galaxies at z^ 2 accrete signifi- 
cant amounts of cold gas, which after accumulating and 
forming gaseous disks, will become gravitationally un- 
stable and lead to the observed high sta r formation rates 
(jGenzel et alJl2008t IDekel et al. Il2009bft . In this picture, 
mergers and hydrodynamic processes like feedback from 
star formation and AGN play only a minor role in assem- 
bling early galaxies, beginning to play a significant role 
only as a consequence of the build-up of the bulk of the 
stellar mass in bulges and disks. 

Given the faintness of the targets and the often low 
number of diagnostic optical emission lines, usually only 
studying Ha or [OIII]A5007, many of these results must 
rely on the assumption that the gas conditions in the 
interstellar medium will overall be largely similar to those 
in galaxies at low redshift. This assumption has not been 
tested directly on the observed properties of high-redshift 
galaxies. 

We already know of (or may plausibly expect) several 
major differences between galaxies at high and at low 
redshift, which may strongly influence the state of their 
interstellar medium and thus, their rest-frame optical 
line emission. First, many galaxies at high redshift are 
less evolved, with hig her fractions fErb et al. Il2006ch of 
lower - metallicity gas (|Erb et al. I l2006at iMaiolino et al.l 
2008). This may have rather subtle observational c onse- 
quences. For example, iRobertson fc Bullockl (|2008f ) dis- 
cuss the expected morphology of an advanced (or mostly 
relaxed) merger of two gas-rich disk galaxies as would 
be observed using state-of-the-art Integral-Field Units 
(IFUs). They find that their model reproduces the ob- 
served properties of high redshift galaxies as well as pure 
(and isolated) disk models, which are often favored by 
observers in contrast to on-going or advanced mergers 
(INesvadba et alj2006aHGenzel et al. l2006tlGenzel et all 
2008: iWrierht, et, alll2008f )~ 

Second, virtually all high-redshift galaxies studied in 
detail with IFUs have very high surface brightnesses of 
the recombination lines. Given the impact of cosmolog- 
ical surface brightness dimming, which is a strong func- 
tion of redshift (oc (z + l) 4 for bolometric luminosities 
and spectral lines, cx (z + l) 5 for broadband photome- 
try due t o the addition al 'stretching' of the continuum; 
see, e.g.. iTolmanl [19 30) . and current observational lim- 
its, all observations at high redshift will naturally be bi- 
ased towards the most luminous emission-line regions, 
and the most luminous galaxies. As a result, all galax- 
ies so far studied with IFUs have star-formation rates 
of several 10s of M© yr" 1 or more, and typical star- 
formation intensities, SFI, well abo ve the critical t hresh- 
old of SFI cri4 ~ 0.1 M Q yr" 1 kpc~ 2 (|Heckmanll2003D nec- 
essary to drive vigorous outflows in local galaxies. Galax- 
ies with SFI J> SFI cr it are observed to create and main- 
tain strongly over-pressurized bubbles of hot gas from 
the thermalized ejecta of supernovae and massive young 
stars, which expand perpendicula r to the disk plane and 
produce galactic-scale outfl ows ([Heckman et all 119901 : 
iLehnert fc Heckman I Il996allbl ). These outflows may 
play an important role in rendering a starburst "self- 



regulating" (|Dopita fc Rvder lll9^ISiik1ll99l . The ob- 
servational signatures of starburst-driven winds, such as 
characteristic velocity offsets between the r est-frame UV 
absorption lines and th e system ic velocity (|Pettini et al.l 
[20001: lErb et al. I I200H l2006al) . or blue wings in rest- 
frame optical emission lines (e.g. iNesvadba et al 1 20071 
are commonly observed at high redshift. In a detailed 
study of a massive, z~ 2.6 starburst and the related 
outflow of a submillim eter-selected galaxy at z^ 2.6, 
INesvadba et al. I (|2007l) found that overall, the physical 
properties of maximal starbursts at z~2— 3 appear very 
similar to local starbursts. The starburst appears to be 
self-regulating, and characteristic, density-sensitive line 
ratios suggest that the pressures in the starburst region, 
which are ultimately driving the observed outflow, are 
very similar to those in low-redshift starbursts. 

This may have a non-negligible impact on our in- 
terpretation of the physical processes in high-redshift 
galaxies, because much of our well-established low- 
redshift emission line diagnostics relies on the physi- 
cal gas conditions in rather subtle ways. For example, 
when relating the ratios of strong, low-ionization neb- 
ular emission lines with Ha, like, [Nil] /Ha, [Oi]/Ha, 
or [Sn]/Ha with the [Om]/H/3 rati os, starburst galaxies 
will fall onto a characteristic curve ([Baldwin et al.|[l981l : 
iVeilleux fc OsterbrocU[l987l ). whereas AGN will fall into 
a different part of the diagram. For high-redshift galax- 
ies, ho wever, the s ame re lationship may not be strictly 
valid. lErb et al. I (|2006afl observed that, although UV 
selected galaxies at z~ 2 follow the overall shape of 
the starburst-curve in the [Nil] /Ha versus [Om]/H/3 di- 
agram, as_j 1 _wdiol_e J _they are_offset towards larger [Nil] 
fluxes. iBrinchmann et a l. (2008) collected a comparison 
sample of low-redshift SDSS galaxies with similar off- 
sets, and found that these galaxies also had statistically 
higher Ha equivalent widths. They argue that the most 
likely explanation may be higher ionization parameters 
and densities in high-redshift star-forming regions com- 
pared to local galaxies. We will further develop these 
arguments and illustrate that pressures induced by the 
starbursts may very naturally explain some of the kine- 
matic properties of high-redshift galaxies as well. 

However, the line ratios indicative of excitation by a n 
AGN may also shift at high redshift (|Groves et al.ll2006h . 
Since most luminous AGN at l ow redshift preferentiall y 
reside in massive galaxies (e.g.. iKauffman n et al 1 120031 
and because gas metallicit y correlates with galaxy mass 
(e.g.. lTremonti et alll2004l ). most AGN hosts at low red- 
shift wfiJ_Jiave_jMrrowdine regions with high metallici- 
ties. iGroves et all (|2006l ) modeled the diagnostic line ra- 
tios for AGN in relatively low- metallicity host galaxies. 
They showed that star-formation plus AGN in relatively 
low metallicity host g alaxies could explai n the offset ob- 
served for example by lErb et aT~l (|2006af ) . We will show 
that many z^ 2 — 3 galaxies fall very close to these re- 
gions, which may make it difficult to robustly quantify 
the role of AGN and starbursts in exciting the optical 
emission line gas. 

We will in the following present an analysis of 11 
galaxies with particularly deep near-infrared integral- 
field spectroscopy obtained at the V1T. These data allow 
in particular to trace intrinsically fainter nebular emis- 
sion lines such as [Oi]A6300 and [Sn] AA6716,6731 which 
being close in wavelength to Ha often fall into the same 
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Fig. 1. — The various emission line and kinematic maps for the source, ZC782941, used in our analysis, (top row and from left to right) 
The Ha flux, velocity dispersion, velocity, [NIIJA6583, and [Nil] A6583/Ha ratio maps respectively for the 250 mas (non-AO assisted) data 
for ZC782941. (bottom row) The maps of ZC782941 in the same order as for the top row but now for the observations taken with the 
100 mas pixel scale and with the assistance of adaptive optics. In each of the panels, we show the approximately size of the full-width at 
half-maximum of the seeing disk (the disk marked "seeing" in each plot) and the relative positions are in arc seconds. The contours in 
each maps are the distribution of the Ho emission in each of the two cubes (with the lowest contour and step between contours of 1.4 X 
10~ 16 erg s _1 cm~ 2 arcsec -2 for the non-AO assisted data and 1.7 X 10~ 16 erg s — 1 cm -2 arcsec - 2 for the AO assisted data). All maps 
were centered on the peak in the Ho surface brightness distribution which defines the (0,0) in each panel. 



band. A subset of 5 galaxies also have measurements 
of the [Om]A4959,5007 and H/3 lines, which fall into the 
near-infrared H band, and which have not been discussed 
previously in the literature. Many of the arguments in 
previous papers are based on the implicit assumption 
that most of the emission originates from "ordinary" Hn 
regions similar to low-redshift disk galaxies and thus, 
that observations of the warm ionized gas are represen- 
tative of all of the kinematics and other properties of 
the phases of the ISM. In contrast, our analysis is aimed 
at quantifying the physical conditions in the interstellar 
medium of these galaxies principally by investigating the 
surface brightnesses of the recombination lines, various 
optical emission line ratios, and line widths to test this 
underlying assumption. 

Throughout the paper we adopt a flat Ho =70 km s _1 
Mpc -3 concordance cosmology with tt\ = 0.7 and Hm = 
0.3. 

2. OBSERVATIONS AND DATA REDUCTION 

For our analysis of the physical conditions in 
the interstellar medium of strongly star-forming, 
high-redshift galaxies, we collected a sample of 
11 galaxies at z~1.5— 2.5 fr om the SINS program 
(|Forster Schreiber et al. I l2006f ) with rest-frame optical 
integral- field spectroscopy (Table [l}. 

Data were taken with the near-infrared integral-field 
spectrograph SINFONI on the ESO Very Large Telescope 
in several runs between 2004 and 2006. Observations 
have been presented elsew here (jForster Schreiber et al. I 
120061: iGenzel et al]|2008h . We reduced these data in- 
dependently froni any other previous work, using the 
IRAF (jTodv 1 119931 ) standard tools for the reduction of 
longslit spectra, modified to meet the special require- 
ments of integral-field spectroscopy, and complemented 
by a dedicated set of IDL routines. Data are dark-frame 
subtracted and flat-fielded. The position of each slitlct 
is measured from a set of standard SINFONI calibra- 
tion data, measuring the position of an artificial point 
source. Rectification along the spectral dimension and 
wavelength calibration are done before night sky sub- 
traction to account for some spectral flexure between the 



frames. Curvature is measured and removed using an 
arc lamp, before shifting the spectra to an absolute (vac- 
uum) wavelength scale with reference to the OH lines in 
the data. To account for variations in the night sky emis- 
sion, we normalize the sky frame to the average of the 
object frame separately for each wavelength before sky 
subtraction, correcting for residuals of the background 
subtraction and uncertainties in the flux calibration by 
subsequently subtracting the (empty sky) background 
separately from each wavelength plane. These data re- 
duction procedures have fewer interpolations and are op- 
timized for faint, extended, low surface brightness objects 
compared to the SINFONI pipeline which obviously must 
be capable of reducing a much wider variety of objects. 
Overall, we expect that these differences lead to more 
robust data compared to previous reductions but overall 
the results are largely consistent. 

The three-dimensional data arc then reconstructed and 
spatially aligned using the telescope offsets as recorded 
in the header within the same sequence of dithered ex- 
posures (about one hour of exposure), and by cross- 
correlating the line images from the combined data in 
each sequence, to eliminate relative offsets between dif- 
ferent sequences. Telluric correction is applied to each 
final cube. Flux scales are obtained from standard star 
observations taken every hour at similar position and air 
mass as the source. 

We also used the standard stars to monitor the seeing 
during observations, and we find an effective seeing in the 
combined cubes of typically FWHM 0"5-0"8. The spec- 
tral resolution was measured from night-sky lines and is 
FWHM- 115 and 150 km s" 1 in the K and H-bands, 
respectively, for the 250 mas pixel scale ("mas" is milli- 
arcseconds). In addition, for two of the sources, data 
were taken with adaptive optics assistance which yielded 
a seeing about 200 mas. The data were taken with the 
100 mas pixel scale in SINFONI. In Fig. [1] we show an 
example of various maps that have not been presented 
previously for the source, ZC782941. In addition, we 
show example maps of the sources, Q2343-BX610 and 
Q2343-BX5 28 for comparison with previo usly published 
maps (e.g., iForster Schreiber et al. I [20061 ). In our sub- 



4 



Lehnert et al. 




•01 iq ii ij ■■ -it 40 oi oi ■( -s 



Fig. 2. — The various emission line and kinematic maps for the source, Q2343-BX610 and Q2343-BX528, used in our analysis. We show 
these two sources for comparison with previously published maps in the literature, (top row and from left to right) The Ha flux, velocity 
dispersion, velocity, [NIIJA6583, and [Nil] A6583/Ha ratio maps respectively for the for Q2343-BX610. (bottom row) The maps of Q2343- 
BX528 in the same order as for the top row. In each of the panels, we show the approximately size of the full-width at half-maximum 
of the seeing disk (the disk marked "seeing" in each plot) and the relative positions are in arc seconds. The contours for the maps of 
Q2343-BX610 represent the continuum, while those of Q2343-BX528 represent the Ho surface brightness (with the lowest contour and step 



between contours of 4 X 10 17 erg s 1 cm 
which defines the (0,0) in each panel. 

sequent analysis, we will analyze the integrated H-band 
spectra of the sources when they are available. We show 
an example of the H and K-band integrated spectra for 
the source, Q2343-BX610, in Fig.O 

In addition, we reduced publicly available 
HST/NICMOS images of 5 of the galaxies in our 
sample, which were observed as part of proposal ID 
10924 (P.I. Shapley). Each galaxy was observed for 
4 orbits using the NIC2 camera with the F160W (H- 
band) filter and a pixel scale of 0"075. The calibrated 
individual exposures (_cal files) were downloaded from 
the HST archive and reduced using the standard IRAF 
routines pedsky, to correct the well-studied pedestal 
effect caused by residual bias, and multidrizzle, to 
combine the exposures. Images were corrected for the 
impact of additional cosmic ray events during passage 
through the South Atlantic Anomaly when required. 
The resulting images are shown in Fig. 0) 

3. THE REMARKABLY HIGH Ho SURFACE 
BRIGHTNESSES AND TRENDS 

In this section, we will discuss the emission line proper- 
ties of our sample of galaxies such as surface brightnesses, 
the relationship between the velocity dispersion and sur- 
face brightness, and the impact of "beam smearing" . 
This analysis will form the bas is of our subsequent dis- 
cussio n. We refer the r eader to iForster Schreiber et al. I 
(12001 . lOenzel et all l[20M) . and iGenzel et alJ (l2008h 
for a discussion of the overall kinematic properties of the 
galaxies in our sample. 

3.1. Ha surface brightness 

One of the most remarkable, but so far little discussed, 
observational findings of z^ 2 — 3 galaxies studied with 
IFUs are their high emission line surface brightnesses. 
In our sample we measure Ha surface brightnesses of 
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arcsec -2 (Fig. [5}. The 
data were smoothed by 3x3 pixels, i.e., averaged in ar- 
eas of size 0"375x0"375, which is appropriate since the 
spatial resolution is typically ~ 0"5-0"6 or ~4— 5 pixels. 
Correcting for cosmological surface brightness dimming, 
this corresponds to a typical rest-frame surface bright- 
ness of 10 -132 - 10 -14 - 8 erg cm- 2 s -1 D" (Fig.©. This 



arcsec ). All maps were centered on the peak in the Hq surface brightness distribution 



represents a lower limit to the highest intensity regions, 
due to the "beam smearing" effect of the low spatial reso- 
lution of our data. In fact, for two of the galaxies, we can 
further quantify the effect of beam smearing on both the 
surface brightness distribution and line widths. For BzK- 
15504 and ZC782941, where we have both adaptive optics 
assisted and seeing limited observations, we find that the 
offset in surface brightness between both sets of observa- 
tions is typically more than a factor of 3 (Fig. [7]). This 
factor is similar to the ratio of the area of the point spread 
function in both sets (the full width half maximum of the 
seeing disk is typically about 0"5 for the seeing limited 
data versus about 0"2 in the AO assisted data), suggest- 
ing that the star-forming clumps are at best marginally 
resolved in either data set. Much of the structure must 
be smaller than about 2 kpc, the physical resolution of 
our highest resolution data. Moreover, this implies that 
all of the measured surface brightnesses of the most in- 
tense Ha emission regions are under-estimated, even in 
the AO-assisted observations. 

An Ha surface brightness comparison between galax- 
ies in the local and distant Universe illustrates the ex- 
ceptional nature of the galaxies in our sample. Galaxies 
that reach the highest surface brightnesses we observe 
and over similar physical scales do not exist at low red- 
shift. In a study of 84 Virgo cluster and isolated spi- 
ral galaxies, many galaxies reach surface brightness lev- 
els at the low end of what we have observed at high 
rcdshift (log SEW <^ —14 erg cm -2 arcsec -2 ) but only 
in their nuclei and on scales about or less than 1 kpc 
(IKoopmann et al~l 2006a: Koopmann fc Kennev |[2006bl 
l2004allbl : IKoopmann et al. II2001D . Local starburst galax- 
ies can reach higher surface bright nesses (but not as high 
as our peak surface brightnesses, iLehnert fc Heckman I 
119951 Il996af) but again only in regions that are nu- 
clear or circum-nuclear with sizes <J 1 kpc (and well 
within the "turn-over" rad ius of the rotation curve; 
ILehnert &: Heckman J [l996b) . Even in the more extreme 
starburst sample of lArmus et al. I (|1989l . Il990f l. only a 
handful of galaxies reach surface brightnesses sufficient 
to be observed at z~2 (such as M82). At high red- 
shift, Ha surface brightnesses as extreme as those found 
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Fig. 3. — The integrated spectrum of Q2343-BX610 in the H-band (left) and the K-band (right). All of the strong optical lines are 
indicated. 
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Fig. 4. — HST/NICMOS NIC2 H-band (F160W) images of 5 galaxies within the sample. Images are 5x5", and North is at the top and 
East is on the left. The images have a pixel scale of 30 mas (the drizzled pixel size) and each image reaches a depth of about 28.6-28.7 
uij\B pixel -1 (for a NICMOS pixel scale of 76 mas) or about 25.8-25.9 m^s arcsec - 2 . 



only in the nuclei of nearby galaxies on small scales are 
found over significantly larger isophotal radii of order 
10-20 kpc and are generally more extreme. Although 
our comparison includes some of the most powerful and 
intense starbursts in the local Universe, they generally 
do not reach sufficiently high surface brightnesses over 
large enough areas to correspond to any of the galaxies 
in our distant galaxy sample. Moreover, at the 4—5 kpc 
resolution of our seeing limited high-redshift data, none 
of the local starbursts would reach these maxima in the 
surface brightnesses that are observed due to the heavy 
spatial smoothing and dilution by regions with lower sur- 
face brightness. 

This lack of correspondence renders any simple anal- 
ogy between "ordinary" quiescently star-forming spiral 
galaxies in the local Universe and galaxies at z~2 ques- 
tionable. Possible significant differences include the dis- 



tribution of mass in various phases of the ISM, of den- 
sity, of star-formation intensity, of gas fraction, and many 
more. Moreover, this lack implies that we cannot only 
study the kinematics of z^2 galaxies without investi- 
gating the nature of the physical processes that power 
their high surface brightness line emission, and the im- 
pact such a finding has on their kinematics and overall 
gas physics. 

3.2. Ha surface brightness: Characteristic correlations 
and Beam Smearing 

To elucidate the underlying physical cause of the high 
Ha surface brightness of our galaxies, we searched for 
correlations with other parameters. As shown in Fig. [5l 
[6l and [3 we find that both the Ha surface brightness 
and Ha velocity dispersion decline with radius. We gen- 
erally do not find a substantial increase in line widths 
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Fig. 5. — A plot of the observed Ho surface brightness ver- 
sus the projected physical radius. The surface brightness has not 
been corrected for cosmological dimming. Each point represents 
a position centered on one pixel of 0'.'125x0'.'125 , whereas the 
data were smoothed by 3x3 pixels, i.e., averaged in areas of size 
0'.'375x0'.'375. As the spatial resolution is typically ~ 0'.'5-0'.'6 (or 
~4— 5 pixels), the points are therefore not independent, and each 
dat a set has a limited nu mber of independent points (about 5-10, 
see [Shapiro ct al. 2008). The legend to the figure indicates each 
galaxy in our sample. The zero radius was chosen to represent 
the symmetry point in either the kinematics or the outer isopho- 
tal Ha distribution. The distribution of surface brightness is not 
symmetrical for most of the sources. 

when comparing the seeing-limited data with data taken 
with the adaptive optics system, although as expected, 
we observe higher Ha surface brightness in some regions 
(Fig. [7]). Thus beam smearing does have a tendency 
to lower the observed surface brightness suggesting that 
most of the structure within these galaxies is not re- 
solved. 

However, beam smearing is worrying in that the trends 
we observe will obviously be influenced by our poor res- 
olution. We must take care in determining what the 
true impact of our poor resolution might actually be on 
the distribution and kinematics of the emission line gas. 
We observe the trends with Ha surface brightness over 
physical scales that are ~6x larger than the spatial res- 
olution of our data. To investigate explicitly whether 
this may be an artifact due to low spatial resolution we 
constructed several simple models that have bright cen- 
tral point sources with broad lines. Such models would 
correspond to a bright AGN or concentration of mass 
increasing the velocity dispersion in the centers of these 
galaxies. Although this may appear unphysical, since we 
do not see bright point sources, it will help to investi- 
gate, for example, the effects of a narrow-line AGN. Our 
toy models (not shown) do not reproduce the observed 
trends, which implies it is very difficult to contrive a 
realistic situation where only the poor resolution of our 
data would lead to the trends we observe between surface 
brightness and radius. 

Could beam smearing induce an apparent relationship 
between line width and surface brightness? Could this 
trend be related to distant galaxies having more com- 
plex light profiles? To understand the impact of beam 
smearing and complex light profiles of these galaxies on 



E 



300 


I 


i i i i i i 


i i i + i i 
+ - 


250 






+- 

+ - 
+- 


200 
1 50 


= Z '' ' ' '' ' ! 

- +BlK- 15504 lOOwncn 

4- 

4- 

■ ++ + 

+ + \ * 

+ ++ ^XJ- 


,\ 

&+ + 

"t-4 -H- 


++ w ++ - 
#Er<>. V ' 

Wjj^r i 

)t »L_4- 


100 

50 



■f *+i +4+ \ + ; 


i~* , ... , i 


w \ 



-15.0 -14.5 -14.0 -13.5 -13.0 
log SBho (erg s~' cm"' arcsec~ 2 ) 



Fig. 6. — Plot of the He* surface brightness (now corrected for sur- 
face brightness dimming) versus the observed velocity dispersion. 
See Fig.[5]for details on the 3x3 pixel smoothing that was applied. 
The legend to the figure indicates which symbol represents which 
galaxy in our sample. There appears to be a relationship between 
surface brightness in Ha and the line width for both the ensem- 
ble of sources as well as within individual sources. Because of the 
trend between surface brightness and dispersion, velocity disper- 
sions below the resolution of our data, about 50 km s _1 for the 
sources observed in the K-band (all but one of the sources shown 
here) and about 60 km s~ 1 for BzK-6397 have particularly large 
relative uncertainties. 

our analysis, we constructed simulated data cubes us- 
ing the light distributions observed in NICMOS H-band 
images of five galaxies in our sample (Fig. [4j. We as- 
sumed that the line emission follows the distribution of 
the H-band flux, that the velocity d ispersion is 25 km s" 1 
independent of position or radius (jEpinat et al. I [2009), 
and that the rotation curves and pea k velocities are those 
from lForster Schreiber et al~l (|2006l ) . We did not include 
noise in this display as it results in a scatter plot about 
the assumed constant velocity dispersion, masking the 
trend of some pixels to reach high dispersions (which we 
think is extremely important to make obvious). Thus for 
clarity, we do not show a plot of our analysis with noise. 

As we can see in Fig. assuming typical velocity dis- 
persions seen in local disk galaxies and light distributions 
of distant galaxies and smoothing the m to our resolutions 
does not reproduce the data (see also [Wright et al.ll200cl 
iGenzel et all [20081 lF5rster Schreiber et al. 1 120061) . We 
do see that there are some regions of high dispersion, but 
this makes up a small number of pixels whereas a great 
majority of the regions, regardless of their relative surface 
brightness have dispersions similar to what we initially 
assumed (i.e., 25 km s" 1 ). This, coupled with the lack 
of increase in the dispersions when comparing our seeing 
limited data with that taken using adaptive optics, sug- 
gests that beam smearing, while obviously playing a role 
in these trends, does not cause these trends. 

4. EMISSION-LINE PROPERTIES OF HIGH-REDSHIFT 
GALAXIES 

Given the high sta r-formation rates 

(|Forster Schreiber et all 120061 : IGenzel et all I2008T) 
and emission line surface brightnesses in these galaxies, 
we hypothesize that the high surface brightnesses and 
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relationship between velocity dispersion and surface 
brightness in Ha, may be explained by postulating 
that the intense star formation is pressure-driven by 
mechanical energy input from the starburst itself and 
self-gravity of gas at high surface densities. In this 
sense, the star formation will be self-regulated. We 
will argue that these systems may be analogous to 
nearby starburst galaxies, such as M82, except that 
at high redshift, the intense star formation and strong 
mechanical energy injection must act over significantly 
larger areas (1 0-20 k pc compared to a few kpc, e.g., 
iHeckman etafl 119901 : iLehnert fc Heckman I Il996a| ). but 
with a similar local surface brightness and similarly 
high pressures. However, this is not the only possibility 
and we will explicitly address different scenarios to 
investigate whether these relationship may be generated 
by cosmological gas accretion or gravitationally unstable 
disks. Some of our main arguments will rely on the 
detection of nebular emission lines like [NII]A6583, 
[SII]AA6716,6731 or [OI]A6300, which are relatively faint 
in high-redshift galaxies. We detected all of these lines 
in only one galaxy, while we detect all but [OIJA6300 in 5 
others. Where one of the lines is undetected or severely 
affected by a telluric night sky line, we give upper 
limits, provided that they are physically meaningful. 
For a subsample of 5 targets, our data include H/3 and 
[OIII]AA4959,5007, which fall into the near-infrared 
H-band for redshifts z>2. These data sets have not 
been discussed previously in the literature. We list the 
integrated emission line properties of these galaxies in 
Table 1. As the H-band data are rather shallow relative 
to the K-band, we typically detect only the highest 
surface brightness regions (and likely also regions of 
relatively low extinction). 

We extracted emission line ratios from matched aper- 
tures in the H and K band data sets, covering the areas 
where line emission is detected in the H band. We use 
the measured, and not the extinction corrected line fluxes 
for the emission line diagnostics. This adds only a minor 
uncertainty, given the relatively low signal-to-noise ratio 
of our data, low luminosities of the low-ionization lines, 
and the fact that we will only use ratios of lines with very 
similar rest-frame wavelengths. 

4.1. Extinctions 

We estimated extinctions for the four objects with 
measured Ha and H/3 fluxes (Table 2). H/3 in Q2346- 
BX482, which would have been a fifth galaxy with a 
H/3 flux and extinction estimate, and the emission in the 
north of Q2343-BX389 are severely affected by night sky 
lines (so the measurement is only for the southern part 
of the galaxy), and for three other galaxies we identi- 
fied the areas with reasonably bright H/3 emission, and 
extracted Ha from the same aperture (Table 2). For a 
galactic extinction law, and an intrinsic Balmer ratio of 
F_hq /Fh/3 =2.86, we find extinctions in the range of 
Ay ~ 1—2 magnitudes. This corresponds to correction 
factors of about 2—5 between observed and intrinsic Ha 
luminosities. We will neglect extinction in much of our 
subsequent discussion, but caution that intrinsic values 
are strict lower limits and may have been underestimated 
by factors of a few. 

4.2. Electron densities and pressures 



The line ratio of the [Sn] doublet, 
[Sii]A6716/[Sii]A6731, is density-sensitive in the range 
of ~ 10 1-5 cm~ 3 . These lines are relatively faint, but 
we have robustly detected and spectrally resolved this 
doublet in six galaxies in our sample, at signal-to-noise 
ratios ~ 10. We are thus able to measure electron 
densities directly, and to estimate the pressure in the 
partially ionized zones within the interstellar medium. 
Densities listed in Tableware in the range ~ 100—1000 
cm~ 3 . Such values are typical, if not higher, than 
densities in the starburst r egions of galaxies in the loca l 
Universe (such as M82; iLehnert k. Heckman I [l996a) 
and suggest thermal pressures of about 10 -8 5 — 10 -9 5 
dyne cm -2 , i.e., several orders of magnitude higher than 
in the interstellar medium in normal nearby galaxies. 
Similar pressures are found for the z=2.6 submilli meter 
galaxy SMMJ14011+0252 dNesvadba et al. 1 12007? ). We 
will argue below that the ambient medium of our galax- 
ies is highly pressurized, and show that these estimates 
are also consistent with pressures derived directly from 
photoionization models. 

The high pressures may also explain another re- 
markable feature of our galaxies, namely their overall 
low ratios of low-ionization lines like [Sn] AA6716,6731 
and [Nn]A6583 compared to the high Ha luminosities. 
While the recombination lines increase linearly with in- 
creasing density and ionization parameter, the ratio of 
[Sn] AA6716,6731 will decline with increasing ionization 
parameter (e.g.. iWang et al. 1I1999D . Since the gas with 
the highest surface brightness also has a declining ra- 
tio, measuring the density from the [Sn] lines becomes 
more difficult as the surface brightness increases. Thus, 
it may not be surprising that we obtained [Sn] measure- 
ments at sufficiently high signal-to-noise only for parts of 
our sample. This may also affect other diagnostic line ra- 
tios. Larger samples of high-redshift galaxies with deep 
spectroscopy of a comprehensive set of faint, diagnostic 
lines will help to secure these findings. 

4.3. Diagnostic line ratios 

iBaldwin et all (|1981h and iVeilleux k. Osterbrockl 
([1987T7 advocated the use of characteristic ratios of the 
bright optical nebular emission lines as diagnostics to 
differentiate between ionization due to starbursts and 
AGN. These so-called BPT diagrams relate the strengths 
of lines like [Nn]A6583, [Oi]A6300, or [Sn] AA6716,6731 
with those of the Balmer recombination lines and 
[Oiii]A5007, and gives us the ability to trace the phys- 
ical conditions, namely temperature and ionization 
parameter, in the emission line gas. Due to the different 
ionizing spectra of starbursts and AGN, galaxies will 
fall into characteristic areas of the diagrams when their 
nebular emission is dominated by photoionization from 
young stars, or by an active nucleus. 

We show the results in Fig. [5] and [TU1 together with 
the loci of other galaxies at similarly high redshifts taken 
from the literature. It should be noted that all of these 
diagnostics are to a large degree empirical and devel- 
oped for galaxies at low redshift, and that different evo- 
lutionary stages may influence the emission line diag- 
nostics in a rathe r subtle way. As already stated by 
lErb et al. I (|20Q6aD . many z^ 2 — 3 galaxies are shifted 
towards higher low-ionizat ion line ratios rela t ive to the 
low-redshift relationships. iBrinchmann et ail (|2008| ) ar- 
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Fig. 7. — Plot of the rest- frame Ha surface brightness versus 
the observed velocity dispersion for the 2 galaxies for which we 
have both seeing limited and adaptive optics assisted data sets 
(ZC782941 and BzK-15504). Individual points in the data cubes 
were treated as shown in Fig. \E\ In the case of the data cubes that 
were taken with adaptive optics assistance, the pixel scale is 50 
mas pixel -1 and they were also averaged over 3x3 pixels (Fig. ??. 
There is an increase in the velocity dispersion of some of the regions 
within BzK- 15504 but these are entirely associated with the region 
of and around the AGN , substantiating our claim that BzK-15504 
hosts an AGN (see also lGenzel et al. |j2006D . 

gued by analogy with a subsample of local galaxies from 
the SDSS that this may be a result of higher pressures 
and higher ionization in high redshift galaxies. 

However, this is not the only effect we may expect 
to take place. At low redshift, luminous AGNs reside 
predominantly in relatively massive galaxies which have 
comparably high (;>solar) metallicities. At high redshift, 
this is not necessarily the case. iGroves et all (|200l mod- 
eled the expected position of AGN in low-metallicity 
host galaxies in classical emission line/ionization dia- 
grams (e.g., [Nil] /Ha versus [OIII]/H/3). They find that 
the metallicity-sensitive [Nil] /Ho ratio will be shifted to- 
wards lower values for AGN with low-metallicity narrow 
line regions. Seyfert galaxies and QSOs have ratios of 
[OIII]/H/3 that are higher than those of HII regions and 
star- forming galaxies in such diagrams. Since the ra- 
tio of, for example, [OIII]/H/3 is not strongly affected by 
lower the metallicity, the position of low-metallicity AGN 
in ionization diagrams may lie above the locus occu- 
pied by HII regions. We illustrate this effect by showing 
the position o f the low-redshift, low-metallicity AGN of 
IGroves et all (|2006f ) in the same diagram. Interestingly, 
most optically or UV-select ed starbursts fall very close to 
the region spanned by the IGroves et al.l sample. This is 
not to say that these will be low-metallicity AGN, in par- 
ticular, since the observed strong starbursts will create 
high-ionization, high-pressure environments, which will 
shift the galaxies toward s the AGN part of the diagram 
(|Brinchmann et ail 2008). However, this does imply that 
we may significantly underestimate the fraction of high- 
redshift galaxies with AGN of rather moderate luminos- 
ity. In the following we discuss the specific example of 
an obscured quasar with relatively inauspicious diagnos- 
tic line ratios compared to optically selected QSOs. 



Fig. 8. — A plot of surface brightness versus velocity disper- 
sion for a set of toy models. Each of the "models" was gener- 
ated assuming the distribution of light of the HST/NICMOS im- 
age, an intrinsic and uniform velocity dispersion of 25 km s — x , a 
distribution of rotation velocit y and a velocity peak as given in 
Forstcr Schrcibcr et al. (2006), and a point spread function with 
a FWHM of 0'.'6. Care was taken to ensure that the final resolution 
of the data set was 0'.'6 and took into account the intrinsic reso- 
lution of NICMOS. For the presentation here, we did not include 
noise as it would simply tend to increase the scatter, especially at 
the low surface brightnesses, but not change the overall trend that 
the dispersion would remain roughly constant. 

5. THE POWERFUL ACTIVE GALACTIC NUCLEUS OF 
BZK-15504 

We will now give an example from among the galax- 
ies of our sample, BzK-15504, where the AGN plays 
a role that is almost certainly non-negligible in inter- 
preting the extended emission of the host galaxy. This 
galaxy falls within the range of low-metallicity AGN in 
the [Nil] / Ha versus [O l ll]/H 3 diagnostic diagram dis- 
cussed by IGroves et al.l (|2006l ) , but also near the divid- 
ing line between AGN and starbursts. Thus, the im- 
pact of AGN photoionization will easily be missed with 
a diagnostic based only on local galaxies since power- 
ful AGN embe dded in low metalli city host galaxies are 
relatively rare ([Groves et~alll2u06l) . We use additional 
constraints, in particular the [Oi]/Ha ratio and the near- 
nuclear [Oiii]A5007 luminosity to argue that this galaxy 
hosts an obscured quasar, and that most of the extended 
Ha emission may in fact be part of an AGN ionization 
cone. This calls into question previous interpretations of 
the star-formation properties and accretion rates in this 
galaxy (see also § 18. 2p . While this is the only clear case 
among the 11 galaxies discussed here, it does illustrate 
that the rest-frame optical emission line diagnostics for 
z^ 2 galaxies, including the 11 discussed here, may be 
less straight-forward to interpret than previously real- 
ized. 

5.1. The Bolometric Luminosity of the QSO in 
BzK-15504 

It is clear that the nuclear emission lines of BzK-15504 
lie within the region of the AGN in all of the emission line 
diagrams suggesting that it is indeed a powerful AGN 
(Fig. [5] and ITU)) . The bolometric luminosity is an im- 
portant parameter as it sets the total energy output of 
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Fig. 9. — The optical emission line ratios, [Nil] A6583/Hct ver- 
sus [Om] A5007/H/3 for some of the galaxies in our sample (BzK- 
15504, BX610, BX482, MD41, and BX389), as w ell as for some dis- 
tant and local galaxies taken from the literature IIErb et al. 12006a ; 
Groves et al.l2006tlNesvadba et al. I2007l;lvan Dokkum et al.l200E ; 
Takat a et al. 1 120061) . The dividing lines within the diagram de- 
marcate the separate regions occupied by galaxies whose emis- 
sion lines are dominated by heating due to star formation from 
those dominated by active galactic nuc lei jKauffmann ct al~l l2003l ; 
Kewlcy ct al. 2001; Osterbrock 1989). Interestingly, many of the 
galaxies in this sample fall near the dividing line, which is consis- 
tent with either heating due to AGN, as is the case for both the 
nucleus and the extended emission line gas in BzK-15506. 

the AGN. BzK-15504 has not been observed over a suffi- 
ciently wide range of wavelengt hs to estimate it s bolo - 
metric luminosity accurately. iHeckman et al. I (|2004f ) 
argue that [Oiii]A5007 luminosity can be used to esti- 
mate the bolometric luminosity of AGN over the range 
of L[ O/// ]=10 6 5 to 10 9 L . Their adopted relationship 
is Lboi =3500L[Q HI ] wit h a scatter of about a factor of 
2 (see IHeckman et al. 1 12004 and references therein for 
details). Using this relation for our [Oiii]A5007 flux from 
BzK-15504 centered on the brightest continuum and line 
emitting region suggests that it has L[ O/// ]=10 9 ' 5 L Q and 
an implied Lb o z=10 13 L Q . Extinction correcting the 
[OlIl]A5007 flux would increase it by 0.5 dex. BzK-15504 
hosts a powerful AGN indeed - a QSO! 

5.2. BzK-15504 as a giant Narrow Line Region? 

In addition to the nuclear line ratios, the line ratios 
of the extended emission in BzK-15504 suggest that it 
could be photoionized by the AGN. In Fig. [H] and [TU1 
we show the BPT diagrams for various line emitting re- 
gions within BzK-15504, and the line ratios lie either near 
the star formation- AGN boundary, as in the diagrams of 
[Olll]/H/3 versus [Sli]/Ha and [Olll]/H/3 versus [Nil] /Ha, 
or clearly within the AGN region as in the [Om] /H/3 ver- 
sus [Ol]/Ha diagram. 

This of course is perhaps obvious from the fact that 
the surface brightness of the extended emission line re- 
gion of BzK-15504 is very high, one of the highest in the 
sample (Fig. [5]). However, can illumination by a central 
AGN explain the light profile of BzK-15504? The light 
profile of BzK-15504 has a ^ in its Ha surface brightness 
dependence and is thus consistent with photoionization 
from a point source such as an AGN. However, since the 
emission is somewhat complex, this is not a strong con- 
straint as other light profiles might equally well fit the 
data. Our point here is to suggest that it is at least "con- 
sistent" with a 4j profile. Making this assumption allows 
us to make a rough estimate of the ionization parameter, 



Uo = S Q so/cn H ; n ff = (1) 

where c is the speed of light, Sqso is the photon in- 
tensity at radius r from the QSO, and where we used, 
Lqso = 10 45 5 erg s" 1 (or 1/10 of the total bolometric 
luminosity), a radius of about 6 kpc, a rough density of 
500 cm -3 and a factor of 3 between the electron den- 
sity and total density to reflect the fact that most of the 
[Sn] AA6716,6731 forms in the partially ionized zone. Us- 
ing these numbers suggests that U is about 10~ 3 . As we 
will see later, this is a rather canonical number for some 
of the other galaxies in the sample and not surprising 
given the surface brightness of all the galaxies. 

Looking at this from another perspective, we can es- 
timate that it takes a few percent of the total bolomet- 
ric luminosity to explain the total Ha emission in BzK- 
15504, if we make a simple recombination estimate for 
the number of ionizing photons. This is less than the 
10% we assumed above, but overall consistent with the 
nebula being completely powered by the QSO. We note 
that this is obviously a lower limit as we have not con- 
sidered the extinction in the extended nebular emission. 

The total ionized gas mass is straight-forward to esti- 
mate as well. If we assume simple case B recombination 
as we did previously to estimate the total ionizing lumi- 
nosity necessary to power the nebula, we find that we 
need about, 

M H = 7n m v n e 

hv Ha a£ J a 

= 9.73 x 10 8 i HQ ,4 3 n-^ 00 Mq (2) 

Using ?i e =500 cm" 3 and I//f cti 43=10 ' 8 , we find a total 
mass necessary to sustain the Ha luminosity of about 
3 x 10 8 Mq. This is a rather modest amount of mass 
and shows that it is very simple to have the AGN ionize 
this much gas, which is only a small fraction of the total 
mass. 

Given these estimates, we appear to be able to ex- 
plain the emission line ratios with our estimated ioniza- 
tion parameter, looking at the emission line diagrams 
and physical parameters for photoionized nebulae in 
I Groves et all (|2004 I2OO60 . The line ratios are indeed 
consistent with having Hydrogen number densities of- 
order_ 100-1000 cm" 3 and a dilute ionizing energy field 
(log Uo ~ -3). 

There are of course analogs to this situation in both 
the local and high redshift Universe. At high redshifts, 
optical emission line ionization cones over scales of kpc 
to 100 kpc are seen in powe rful radio galaxies (e.g., 
iNesvadba et al. Il2006bl l2008bl ). Emission line nebulae 
this large are generally seen in AGN with UV and emis- 
sion line luminosities that are higher than observed for 
BzK-15504 (about a factor of a few to 10). However, the 
extended emission line region in BzK-15504 is consistent 
with luminosities observed in the " narrow line regions" 
of distant QSOs (jNetzer et al.ll2004f ). 

AGN also show strong asymmetries in their emission 
line distributions. This may also explain the asym- 
metries seen in the emission line images of BzK-15504 
(|Genzel et al. I l2006f) . It represents the asymmetry in 
the opening angle out of which the photons escape the 
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AGN and where the gas falls within the beam. A local 
example of this is NGC 1068 (jVeilleux et al.l f2003h . In 
NGC 1068, the circum-nuclear ionization cone shows a 
strong asymmetry in its distributio n of [Olil]A5007 on 
both small scales (100s of parsecs; Evans e t al.lfl993h 
and large scale (kpc scales; IVeilleux et al.ll2003l ). Only 
in the areas of intense Ha emission are the line ratios 
consistent with photoionization of massive stars in Hn 
regions. Over a much larger scale, 10s of kpc, in the 
areas of relatively low Hydrogen recombination line sur- 
face brightness are the lin e ratios consistent w ith pho- 
toionization by the AGN (|Veilleux et alj 120031) . More- 
over, the emission line widths over the regions excited by 
the AGN (FWHM~100-few 100 km s" 1 ) are also in rea- 
sonable agreement with what is observed for BzK-15504. 
In many optically selected AGN, the line widths are not 
strongly influenced by the AG N despite the fact that th e 
ionization state of the gas is (|Nelson fe Whittle lfl996h . 
Of course, both the gas mass and the power of the AGN 
are larger in BzK- 15504 and it may well be that as a 
result, the AGN can even outshine any star formation in 
the extended optical emission line emitting gas (Fig. [9] 
and [TO]) . Within this regard it is also important to note 
that BzK-15504 also has one of the highest Ha surface 
brightnesses in our sample. The AGN could be respon- 
sible for increasing its overall surface brightness as well 
as its AGN-like extended line emission. And like many 
classes of AGN (and in analogy with NGC1068 as just 
discussed) the radiation, although likely heavily atten- 
uated, is able to escape to large distances despite the 
ionized gas representing a small fraction of the total gas 
mass of BzK-15504. 

This last point is important. AGN host galaxies show 
a reasonable correlation between the si ze of their nar- 
row l ine regio n and the luminosit y in H/3 ([Bennert et al. I 
120021 but see. lNetzer et al.ll2004[ ). The H/3 luminosity of 
BzK-15504 would suggest a narrow line region size of 
a few to 10 kpc (this is without extinction correction). 
Interestingly, this is approximately the size of the emis- 
sion line nebulae observed in Ha and [OIII]. Thus, in 
agreement with other arguments, there is evidence that 
the extended emission line region is nothing more than 
a narrow line region around a powerful QSO with line 
widths that are, by analogy with local AGN, influenced, 
perhaps dominated, by other processes. 

6. NATURE OF THE EMISSION LINE GAS IN THESE 
GALAXIES 

The surface brightness in our sample is both very high, 
and a function of radius, velocity dispersion, and low 
ionization line emission ratio (i.e., [Nil] /Ha). This is 
in sharp contrast to what is observed for spiral galax- 
ies in the local Universe where the velocity dispersion is 
roughly constant as the surface brightness declines expo- 
nentially. With the local trends in mind, we hypothesize 
that the trends seen in the distant galaxy data can be 
explained by self-regulated star formation whereby the 
intense star formation is pressure driven by the gas mo- 
tions induced by the star formation itself. We attempt to 
show that essentially, these systems are analogous to the 
nearby starburst galaxies, such as M82 (and regulated 
perhaps in the same way as the ISM of the Milky Way) 
with an overall similar surface brightness and similarly 
high pressure, but with their star formation occurring 



over a much larger area. 

6.1. The relation between Ha surface brightness and 
[Nil] /Ha 

In Fig. QTJ we show the relationship between the emis- 
sion line ratio [Nil] /Ha and the surface brightness of Ha. 
An obvious way to understand such a relationship, in ad- 
dition to the high Ha intensity is through the gas pres- 
sure in the cloud or cloud interfaces that gives rise to the 
recombination line emission. 

To test this hypothesis we ran some simple photoion- 
ization models using the code Cloudy 3 . The input spec- 
trum was generated using Starburst99 (jLeitherer et al. 1 
1999) for a constant star formation rate and an age of 
10 s years (consistent with t he estimated ages of galax- 
ies at z^2 lErb et al. I l2006bh However, the exact age is 
not very important as long as it is sufficiently long for 
the spectral energy distribution, especially in the UV, 
to reach an equilibrium shape. It is thus appropriate 
for ages greater than about a few 10s Myrs and con- 
stant star-formation rate. We assumed a constant den- 
sity slab with ionization parameters ranging from about 
log U=— 5 to —1, initial densities log n=l, 2, and 3, and 
ISM metallicities and grain depletion which are kept con- 
stant for all calculations. These attempts at modeling 
are not intended to be exhaustive but illustrative. To 
gauge the impact of the ionization of the modeled cloud, 
we ran these calculations for two total column densities, 
10 20 and 10 21 cm . The results from this modeling are 
shown together with the data in Fig. QT] 

We can see that the range of [Nil] /Ha and surface 
brightness can be explained by a combination of column 
density, volume densities, and moderately diffuse radia- 
tion fields. Obviously such simple assumptions are not 
going to "fit" the data in any sense of the word, but 
show that the area of line ratio-surface brightness space 
covered by the data can be explained by a combination 
of these parameters, with each galaxy exhibiting some 
range in each of these. Surface brightness is linearly pro- 
portional to the density and ionization parameter, while 
the line ratio depends on roughly one over the square 
root of the ionization parameter for the range of ioniza- 
tion necessary to explain the line ratios. Such a high 
density and ionization rate implies that the ISM of these 
galaxies must in general exhibit high thermal pressures, 
P/k~10 6-7 K cm~ 3 or more. The true maximum pres- 
sures are likely to be higher as our physical resolution is 
only about a few kpc making it difficult to identify the 
regions of highest surface brightness. 

These pressures are much higher than observed in the 
disk of our Milk y Way or other nearby norm al galaxies. 
As discussed in iBlitz fc Rosolowskv I (|2006l ). and refer- 
ences therein, the hydrostatic mid-plane pressure in the 
MW is about 10 3 - 3-4 - 3 cm~ 3 K, while for other local spi- 
rals it can range up to about 10 6 cm~ 3 K but is typically 
about 10 46 cm -3 K. There have been fewer estimates 
of the more appropriate comparison, namely the ther- 
mal pressures. In the nuclear regions of local starburst 
galaxies, the thermal pressures as estimated from the 
ratio [Sn] A6716/[Sn]A6731 give values that range from 
10 -9 ' 3 to 10~ 85 dyne cm -2 or typically about 10 7 cm -3 

3 Calculations were performed with version 07.02.02 of Cloudy, 
last described by Fcrland ct al. (1998). 
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Fig. 10.— The optical emission line ratios, [Sll] AA6716,6731/Ha versus [Om] A5007/H/3 (left), and [Ol]A6300/H« versus [Om] A5007/H/3 
( right ) for the same subsample as in Fig. [9] Again, the lines demarcate the regio ns of the emission line diagrams where the photoionization 
is dominated by star formation or AGN( Kewley et al. 2001: Osterbrock 1989). These diagrams again demonstrate the role of the AGN 
in photoionizing BzK-15504. The other galaxies in the sample lie close to the lines of demarcation in the [Sll]AA6716,6731/Hct versus 
[Om] A5007/H/3 diagram, but have limits that are consistent with being photoionized mainly by massive stars. 



K (|Lehnert fc H cckmanl ll996a[) similar to that o b served 
here at z~2. Interestingly, IBrinchmann etTai] ((20081) 
make a similar argument for galaxies that generally fall 
off the locus of the Hn regions in the [Om]/H/3 versus 
[Nil] /Ha diagram and lie close to the dividing line be- 
tween Hn region and AGN excitation. They find that 
such galaxies have high specific star-formation rates, as 
do the general population of high redshift galaxies dis- 
cussed here, and r elatively large Ha equivalent widths. 
iBrinchmann et aTl suggest that the most likely explana- 
tion for this is higher density (pressure) and that the es- 
cape fraction may be higher in the star-forming regions 
implying that the nebulae are at least partially density 
bounded. Similarly, in Fig.QTJwe see that the clouds with 
column densities of log N# =20 become density bounded 
at high ionization parameters and may explain the sur- 
face brightness limits we observe. Thus the pressures 
estimated from the photoionization models are consis- 
tent with those of nearby actively star-fo rming galaxies 
(|Heckman et al.lll990t fWang et al. II1999D . 

6.2. Consistency with high ionization lines 

In addition, for some of the galaxies we have H-band 
data cubes of sufficient signal-to-noise to investigate the 
ratio of [Oiii]A5007 and H/3. The photoionization models 
would predict [Om]A5007/H/3 of about 0.5-2.0. Indeed, 
not considering BzK-15504, which is powered by a QSO, 
the other galaxies show ratios consistent with the high 
ionization line ratios and again amplifying the idea that 
these nebula have incre ased density compared to nearby 
galaxies as observed bv IBrinchmann et al.l (|2008l ). 



7. 



POWERING THE LOCAL MOTIONS IN THESE 
GALAXIES 



We have argued that the ISM of these distant galax- 
ies is under high pressure and that this high pressure 
could be induced by the intense star formation within 
these systems. Fundamentally, this argument is analo- 
gous to the situation in the most intense starbursts in the 



local Universe such as M82 (e.g.. iHeckman et al.l fl990t 
iLehnert fe Heckman |[l996at iLehnert et al. II1999F . Such 
a hypothesis nicely explains some of the unusual features 
in these distant star-forming galaxies such as their high 
surface brightnesses in the recombination line(s) (which 
makes them observable in the first place), the low ioniza- 
tion line ratios, the ratio of the [Sn] lines (suggestive of 
high densities), and the trend for more intense star for- 
mation to lead to broader lines. We observe low ratios 
(r^2-4) of v/er in the extended line emitting regions of 
our objects, much lower than the ratio of ^10 generally 
observed in local and intermediate redshift disk galaxies, 
and this appears to be driven mainly by unusually high 
velocity dispersions, not low velocity shears. 

In principle, the high velocity dispersions observed in 
the ionized gas lines are likely not representative of the 
turbulence in the whole star-forming ISM. Nevertheless, 
the energy injection necessary to sustain the high veloc- 
ity dispersions observed in the ionized gas should affect 
the atomic and molecular phases (although if it is a tur- 
bulent cascade, with energy injection on large scales, the 
velocity dispersion will likel y be lower in denser gas but 
the p er 2 will be the same; Uoung. Mac Low, fc Brvaiil 
2008). Indeed, high dispersions in the gas out of which 
stars form is suggested by the large sizes and masses of 
high-redshift star -forming regions (giant clumps e qualing 
the Jeans mass, lElmegreen fc Elmegreenll2005l) . This 
is also supported by th e thickness of high-redshift disk s 
when observed edge-on (Elmegre en fc Efmegreen"ll2006l) . 

What is the power source of the high velocity disper- 
sions that are observed, and that should affect the whole 
ISM of high-redshift galaxies? There are several possi- 
ble mechanisms: (1) the conversion of potential energy 
of infalling material in turbulent gas which has suffi- 
cient angular momentum to r elax into a disk configu- 
ration (cosmological accr etion; iForster Schreiber et al. I 
l2006HDekel et al. |[2009bT) ; (2) peculiar motions in unsta- 
ble disks that lead to perturbations in the velocity fields 
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Fig. 11. — The Ha surface brightness (corrected for surface 
brightness dimming) versus the logarithm of the [Nil] A6583/Ha line 
ratio for all the individual spectra of the galaxies in our sample. 
The galaxies are labeled in the legend to the figure. The lines 
represent results for photoionization modeling (see text for details) 
for 6 set of conditions. The lines represent a range of density (log 
uh =1 to 3) and column densities of log Nh=21 cm -2 (solid lines) 
and log Nj^=20 cm -2 (dotted lines). The ionization parameters 
span from log U=— 5 to —1.0 (which increases from left to right 
along the lines, meaning low ionization parameters have relatively 
high log [Nil] /Ha and low Ha surface brightnesses). 

that are not resolved in our data and lead to a high a 
without consisting of real small-scale turbulence (e.g., 
iBournaud et al. I l2008f): (3) self- gravitationally gener- 
ated turbulence (jWada et al.ll2002f k and/or (4) the com- 
bined effect of the intense star formation. 

To characterize the ordcr-of-magnitude needed to 
power the turbulence, we shall start by estimating its 
energy budget (if that is what these motions represent). 
Our low spatial resolution in the rest-frame of the galax- 
ies does not allow us to cleanly separate these possible 
sources of the o bserved broad line emission. In fact, 
iMac Low"! (|1999f ) as well as our astrophysical reason- 
ing below imply that bulk motion and turbulence are 
intricately related in some of the scenarios, in the sense 
that bulk motion will inject kinetic energy into the sys- 
tem which will then be dissipated through turbulence. 
However, carefully examining the total energy contained 
within the emission line gas may help us constrain the 
mechanism responsible for its characteristics. We first 
consider the energy dissipated within the turbulence and 
whether or not it is consistent with cosmological accre- 
tion of gas. 

7.1. Turbulent energy dissipation 

The violent motions observed in these distant galaxies 
are highly supersonic given the densities derived previ- 
ously for the optical emission line gas and the turbulence 
generated, compressible. To estimate the dissipation of 
the turbulent energy therefore requires comparison with 
(magneto-)hydr odynamic simulati ons of a realistic inter- 
stellar medium. IMac Lowl (|1999h provides an estimate 
of the energy dissipation rate for compressible turbulence 
as, 

Efcjn ~ -n u mkvf ms (3) 

where Uy is a constant of proportionality, which 
IMac Lowl ()1999l ) estimates to be 0.21/w, v rms is the root 



mean square of the velocity in the region, m is the to- 
tal mass, Ekin is the kinetic energy, and k is the driving 
wavenumber. Although, the validity of this energy dis- 
sipation estimate on galaxy scales has yet to be verified, 
it is useful to give an order-of-magnitude estimate of the 
energy injection rate necessary to sustain the motions we 
observe. 

The parameters necessary to estimate the energy dissi- 
pation rate of compressible turbulence span a wide range 
of values. For example, the velocity dispersions observed 
in the galaxies range from about <40 km s _1 , in the 
outer regions, to about 250 km s -1 in the inner, circum- 
nuclear regions (Fig. [6|). The typical velocity dispersion 
is of-ordcr ^100-150 km s _1 . The driving scale is par- 
ticularly difficult to estimate since the driving mecha- 
nisms likely operate over a wide range of scales (e.g., 
Uoung. Mac Low. &; Brvaril l2008f ). For example, if the 
largest scale driving mechanism was the cosmological ac- 
cretion of gas, we w ould expect this scale to be large 
(|Dekel et al. I l2009bf ) approximately that of the thick- 
ness of the disk. Alternatively, if the driving mechanism 
is star-formation, the size of the largest star-forming re- 
gions or associations in the galaxies might be the ap- 
propriate scale. The clumps seen in Fig. 0] are approx- 
imately 100s of parsecs to kpc in diameter. Since the 
driving wavenumber is inversely proportional to the driv- 
ing length, assuming a small scale for the driving length 
would tend to increase the energy dissipated through tur- 
bulence. 

We cannot estimate the total mass of gas in the galax- 
ies in a straight forward way. For simplicity, we w ill ap- 
ply the Schmidt-Kennicutt law (|Kennicutt lll998ah to the 
star-formation intensities estimated from the Ha surface 
brightness distribution. This relation implies that the 
gas surface densities are of-order E ffas = 10 2 6 M Q pc -2 
for star-formation intensities of 1 M Q yr _1 kpc~ 2 (1 
Mq yr _1 kpc -2 is the typical average value for the 
star-formation intensity approximately averaged over the 
isophotal radius). This of course assumes that the molec- 
ular gas has same kinematics as the warm ionized gas. 
The high densities that we found in § 16.11 and the course 
angular resolution of our data which averages the kine- 
matics over a large region p erhaps imply that this is not 
a bad assumption (but se e Uoung. Mac Low, fc Brvanl 
[20MlWalter et al. l[200l ). 

If we adopt, v rms =150 km s" 1 , a mass surface den- 
sity of 1000 M Q pc -2 , a driving length of 1 kpc (which 
is a typical thickness of the "c lumpy disks" observed at 
simil ar redshifts as our sample; lElmegreen fc Elmegreen 1 
2006), we find that turbulence likely dissipates about 
10 42 erg s _1 kpc -2 . The surface area within the isophotal 
radius is approximately 200 kpc 2 (Table 1) giving a total 
energy dissipation of about 10 4 erg s" 1 . We empha- 
size given the uncertainties, the limits of our theoretical 
understanding of turbulence and untested assumptions, 
this estimate of the turbulent dissipation should be con- 
sidered as order-of-magnitude only. 

7.2. Cosmological gas accretion 

We can compare this rough energy dissipation estimate 
with that of gas infalling on to the disk itself. The total 
energy accretion rate fr om gas falling onto the d isk is 
given approximately by (|Dekel &: Birnboim Il2008f ). 
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where | A<j> | is the potential energy of infall, M gas l oo 
is the halo gas accretion rate in units of 100 M Q yr _1 and 
V c ,2oo is the circular velocity of the halo in units of 200 
km s _1 . Thus it appears that accretion of gas in itself 
cannot power the turbulent and bulk motions we observe 
in these gal axies if these motions decay as compressable 
turbulence. iDekel et al. I (|2009al ) also suggest that the 
main source of turbulence is not infalling gas, unless this 
infalling gas is itself also highly clumpy, something they 
themselves rule out as highly unlikely for a large fraction 
of the infalling gas. 

7.3. Velocity dispersions in Jeans unstable clumps 

The dynamics of these galaxies could be influenced 
by the mode of star formation. Elmegreen and collab- 
orators have proposed that the large number of clump- 
cluster and chain galaxies observed at high resolution 
in deep HST image ry represent gas- rich Jeans unsta- 
ble disks (review in lElm egreen 1 120071 ). Clumpy galax- 
ies are not rare. Their frequency and relatively fast 
dynamical evolution suggests that perhaps all galaxies 
pass through a "clumpy" stage in their evolution and 
that this process is a natural way of explaining phe- 
nomena like the growth of bulges and nuclear supermas- 
sive bl ack holes (|Bournaud et al. Il2007t lElmegreen et al.1 
I2008al) . HST/NICMOS images of five of our sample 
galaxies (Fig. [4} do show that they have "clumpy" ir- 
regular morphologies consistent with this hypothesis. 

The dumpiness of the disks is hypothesized to be 
driven by Jeans instability, which implies a relation- 
ship between the mass of collaps ing gas and the velocit y 
dispersion within the gas (e.g., lElmegreen "et~alll200l . 
Specifically, the Jeans relationship implies that, 

a gas ~ My 4 GV 2 £^ = 54M^ 9 4 S°^ km s' 1 (5) 

where <j gas is the velocity dispersion of the gas, G is 
the gravitational constant, S sas is the gas surface den- 
sity in Mq pc~ 2 , Mj y g is the Jeans mass in units of 
10 9 Mq, and Xsf.r is the star-formation intensity in 
units of M© yr _1 kpc~ 2 . We have u sed the Schmidt- 
Kennicutt relation (|Kennicutt Ill998bl ) to convert from 
gas surface density to star-formation intensity. We show 
the relationship between the velocity dispersion and star- 
formation intensity for Jeans unstable clumps in Fig. [12] 
for a clump of 10 9 Mq, similar to the largest masses es- 
timated fo£_cbamp_s_b^sed^n_^p^tral_energy distribution 
fitting (Elme green fc Elmegreenl l2005f ). We chose 10 9 
Mq to put the weakest limit on the possible contribution 
of the clumps to the observed velocity dispersion. 

The velocity dispersions predicted by Eq. (6) gener- 
ally lie below the data in Fig. [12l especially at the high- 
est intensities, and the relationship is also too flat as a 
function of star-formation intensity. Although this does 
not rule out clumps in a disk as a contributing source 
to the high velocity dispersions observed in our sample, 
especially for the gas with lowest dispersions, it cannot 
be the entire explanation. And taken at face value, the 



masses necessary to explain the average dispersion would 
be more like 10 10 Mq, not 10 9 M suggesting that the 
dispersions are not dominated by the internal dispersions 
of the clumps themselves. 

7.4. Gravity powering turbulence in dense gas-rich 
galaxies 

The clumpy light distributions observed suggest, in 
particular, that the clumps themselves may have re- 
quired an earlier source of turbulence. This is neces- 
sary to ensure that the initial Jeans mass is as high as 
the clumps masses ob served (up to more than 1O 9 M ; 
lElmegreen et al.l l2007t ) . This implies that we may be 
obliged to hypothesize another source of turbulence, es- 
pecially in initiating the intense star formation, to fa- 
vor and encourage the growth of massive clumps. Per- 
haps self-gravity powering the initial turbulence would 
play the necessa ry role (see iThomasson et al.l 119911 : 
iWada et al.ll200^ . 

A gas disk with a low turbulent speed a may likely 
have a Toomre parameter Q — < 1, if the surface 
density £ is high. As a result, gravitational instabilities 
will form and heat the gaseous medium, increasing its 
turbulent energy until Q reaches ~ 1 and the process 
stops. The observation of clumps in our z~2 galaxies 
and others suggest that Q is close to one, so that this is 
the level at which gravity could indeed power the turbu- 
lence. In local flocculent spiral galaxies. lElmegreen et al.l 
( 2003) argued that gravity was triggering the turbulence 
through local instabilities. 

The hydrodynamic simulations by lAgertz et al. I 
(2009) have shown that disk self-gravity likely triggers 
the 5-10 km s _1 turbulent speed of extended HI disks 
around local spirals - which could not result from star 
for mation beyond the edge of star-forming disks (see 
also IWada et alj I2002D . The density of these modeled 
disks is nevertheless far lower than the estimated den- 
sity of gas disks at z~2. A gas fr action of ~50% in 
these disks (from iDaddi et al. I 20081 observations o r in- 
ferred from the Schmidt law, iBouche et al. I l2007f ) in- 
deed corresponds to typical gas surface densities of ~500 
M pc -2 . More direct hints on the role of self-gravity 
in high-redshift disks c an be found in the mo dels by 
iTasker fc Brronl (|200l ; iTasker fc BrmrTl (f2008h . where 
gas disks with large-scale densities around 100 M Q pc -2 
are modeled, much closer to the hypothesized density at 
z~2 even if still somewhat lower. Interestingly, these 
models show larger clumps when the gas density is large, 
indicating larger Jeans masses and hence larger veloc- 
ity dispersions. This is found in models of isolated disks 
without feedback, where only gravity can power turbu- 
lence. The addition of feedback from star-formation in 
these models does not change significantly the clump 
masses, hence the turbulent speed, suggesting that the 
early stages of turbulence development could be mostly 
powered by self-gravity. Direct models of disks at z~2 
with high density wou ld be desirable, but given that 
ITasker fc Bryan I (|2008l ) have both the gas density and 
the clump masses/sizes somewhat lower than z^2 stan- 
dards, it is reasonable to infer that disk self-gravity is 
likely important source of turbulence at high redshift 
but does not likely generate the extreme dispersions we 
observe (Fig. IT2"]) , As we noted earlier, it may be that 
the formation of massive clumps themselves requires high 
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Fig. 12. — A plot of the star-formation intensity versus the He* 
velocity dispersion in our sample of galaxies. See Fig. 1 for details 
on the 3x3 pixel smoothing that was applied. The symbols for 
the different galaxies in our sample are shown in the legend to the 
figure. The solid black lines show three simple relationship of the 
form a = y/eE, where e is the efficiency of coupling between the en- 
ergy injected and the ISM. (from bottom to top): a = 100y/T,sFR 
km s _1 , <r = IA0^T,sfr km s _1 , and a = 240y/T,sFR km s _1 
(see text for further details). We note that the overall trend seen 
in the diagram for the ensemble of galaxies is also traced by the 
data from individual galaxies but with a somewhat steeper slope. 
If the dispersions represented turbulent motions, we would expect 
a scaling, a = (eE) 1 / 3 . Using two scalings for the coupling effi- 
ciency, 25% and 100% and a primary injection scale of 1 kpc, we 
find cj = SOEg^^j km s -1 and a = lSO^Jp^ km s — 1 (bottom and 
top black dashed curves respectively). The red solid line repre- 
sents the velocity dispersion of a 10 9 Mq clump using the simple 
Jeans relationship between mass and velocity dispersion (see text 
for details). The offset between galaxies may be due to a range 
of coupling efficiencies between the mechanical energy output from 
the on-going star formation, differences in the average pressure 
within the ISM, or geometrical factors such as inclination. 

turbulence initially (Eqn. 5). 

In addition, high velocity dispersions could be a 
result of mergers. In the local Universe, intensely 
star-forming galaxies show very high velocity dis- 
persions in t heir optical emission line gas, up t o 
200 km s" 1 (|Monreal-Ibero. Arribas. fc Colinal 120061) . 
To explain the low v/ a obs er ved at high redshift 
([Forster Schreiber et al.1 l2006h . iRobertson fc Bullock! 
(2008) proposed that a disk settling after a major 
gas-rich merger would have low v/cr. In particular, 
they construct their model to match approximately the 
properties of BzK-15504 in detail including its high 
star-formation rate and intensity (without the AGN). 
Thus, it is not clear how much of the dispersion is 
due to the intense star-formation and how much is 
driven purely by gravity in their mo del. Moreover, re- 
cently |BournauT&^lmegreen] (|2009f ) have criticised the 
merger model in that it does not naturally produce the 
clumpy morphologies that are frequently observed in dis- 
tant galaxies and would lead to disks that make up a 
rather small fraction of the total mass. Thus it is not 
clear if gravitationally driven flows within mergers can 
produce the dispersions and high emission line surface 
brightness we observe without the accompanying intense 
star formation. 



7.5. Ha surface brightness-velocity dispersion: 
Powering the kinematics through star formation 

We argued above that self-gravity may be an important 
source of turbulence at high redshift, but does not seem 
sufficient to generate the observed velocity dispersions 
(see also Fig. IT2|). This and the relationship between Ha 
surface brightness and velocity dispersion suggests that 
the star formation within the galaxies is powering the 
dynamics of the emission line gas. If the star formation 
is indeed inducing the high pressures, then this is what 
would be expected. On purely dimensional grounds, if 
the energy output from young stars is controlling the dy- 
namics of the emission line gas, we would expect, that 
a, the velocity dispersion, would be proportional to the 
square root of the energy injection rate, dE5i?/dt, due 
to stars. If the coupling efficiency of the mechanical en- 
ergy output of the star formation does not depend on 
radius, then the energy injection rate is simply propor- 
tional to the star-formation rate. In this case of course, 
the energy injection per unit area is proportional to the 
star-formation intensity. This hypothesis is equivalent 
to conserving (with some efficiency) the mechanical en- 
ergy output of the star-formation within the ISM of the 
galaxy and that the velocities of the warm ionized gas 
trace this energy injection rate. 

In Fig. we have over-plotted just such a relation- 
ship. This is not a fit to the data, but a simple scaling law 
based on the star-formation rate per unit area and the 
velocity dispersion in the warm neutral/ionized gas in the 
disk of the MW and other nearby galaxies. The function 
is of the form, a=(eSsi?fi) 1 ^ 2 where e has been deter- 
mined for the MW and other nearby galaxies where the 
velocity dispersion in Ha and star- fo rmation intensities 
have been related (see lDib et al.ll2006l . for details). There 
are several possible values for e in this relationship. If 
we take a simple scaling from galaxies like the MW, they 
typically have star-formation intensities in the regime of 
£sftf=10~ 5 to 10~ 3 M yr -1 kpc -2 and velocity dis- 
persions in the warm ionized gas of-order 10 km s . 
iDib et afl suggested that galaxies may change from a qui- 
escent disk mode to a starburst mode at T,sfr=W^ 2 ' 5 
to 10~ 2 M yr" 1 kpc" 2 . This comes from modeling the 
ISM with a coupling efficiency of 25% to the supernova 
energy from disk star formation. Using these two values 
for the scaling relation yields the bottom two curves in 
Fig- [IH whereas a coupling efficiency of 100% would yield 
the curve at the top in Fig. [TSJ We emphasize that the 
axis of abscissa is a simple scaling between Ha surface 
brightness and the star-formation rate which was made 
assuming a rela tionship between st ar-formation rate and 
Ha luminosity (jKennicutt 1 1 1998ah . We have not taken 
into account the effect of extinction in making this scal- 
ing. If the galaxies for which we have H/3 measurements 
for are representative of our entire sample. This would 
increase the star-formation intensities by a factor of a 
few (Table 2). 

While such a toy model is not a particularly good fit, 
especially at the highest star-formation intensities, it has 
the virtue of having no free parameters. It is just a sim- 
ple scaling based on the MW and other nearby galaxies. 
Al though we put t his in context of the modeling done 
bv lDib et all |2006l) . we could as well have simply used 
the MW scaling of star-formation intensity and velocity 
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dispersion. Obviously, the true nature of the interstel- 
lar media of these distant galaxies cannot be so simple 
as the energy injection argument presented here. This 
argument simply shows that the star formation in the 
galaxies themselves is likely controlling the dynamics of 
the emission line gas. 

7.5.1. Mechanical energy due to star formation 

Given the relationship between star formation inten- 
sity and the velocity dispersion of the gas, it is logical 
to investigate whether the star formation has sufficient 
global energy injection rates to drive the high velocity 
dispersions - much like our previous estimate in § 17.11 
comparing the accretion energy rate from i nfalling cos- 
mological gas. Using the relationship from lKennicutt~l 
( 1998h|) between Ha luminosity and star-formation rate, 
we find that the total star-formation rate per unit area - 
the star-formation intensity - ranges from about 0.05-5 
M Q yr _1 kpc -2 (Fig. [5]). The dynamical time and in- 
tegrated Ha equiv alent width of our sample of galaxies 
and that from the lErb et aI71 (|2006bf ) sample, suggests 
that they may have been forming stars continuously at 
the observed rate for the last few x 10 8 years. This 
implies that the ongoing star formation in this sample 
in total produces more than ip 41 - 42 - 5 e rg s _1 kpc -2 of 
mechanical energy (Lcithcr er et al. II1999T ) for the range 
of large scale star-formation intensities observed. These 
numbers are lower limits since we do not have sufficient 
data to constrain the spatially resolved extinction and so 
the true star-formation intensities are likely to be higher 
(Table 2). 

We can now compare this total energy injection rate 
with the value estimated § 17.11 namely 10 443 erg s _1 . 
For a fiducial value of 1 M Q yr _1 kpc" 2 , which is typical 
within a isophotal radius (Fig. [6|) , and the mechanical en- 
ergy output rate given above, would suggest total energy 
injection rates of >10 44 2 erg s _1 . Since the numbers 
are similar, we suggest that the star-formation is power- 
ful enough to maintain the turbulent and bulk motions 
observed. 

7.5.2. Accelerating the emission line clouds 

Having now argued that the star formation is inject- 
ing energy into the ISM of these galaxies and is suf- 
ficient to explain the overall dynamics of the gas and 
to support dissipation through turbulence, can we make 
a plau sible model for how the clouds might be acceler - 
ated? Klein et al. 1 (|1994f ) and INakamura et al. I (|2006l ) 
have developed both analytic models and simulations for 
clouds accelerated in a blast wave. Such a blast wave is 
expected to be generated by the intense star formation 
observed in our sample of galaxies. As the blast wave 
passes it shock heats the cloud and eventually destroys 
it through Kelvin-Hclmholtz and Rayleigh- Taylor insta- 
bilities. 

In this scenario, the turbulence is driven by large 
scale bulk motions induced by energy injection - blast 
waves generated by intense star-formation and other pro- 
cesses - over a wide range of scales, which then cas- 
cades into the denser gas and redistributes the energy 
over all scales in the medium finally dissipating mainly 
on the smallest scales (|j"oung. Mac Low, fc Bryan"ll2008t 
iMac Low fc Klessen Il2004 and references therein) . The 
nature of turbulence emphasizes the importance of both 



bulk and random motions and estimating energy dissipa- 
tion rates based on turbulence arguments while estimat- 
ing velocities based on bulk motion arguments are thus 
a ppropriate. 

iKlein et al. I (|1994l ) have found relatively simple an- 
alytic approximations that can capture the destruction 
time of clouds with a range of properties like contrast 
with the inter-cloud material, Mach number of the blast 
wave, etc. Although there is a problem with simply us- 
ing these formulae with the parameters from the pho- 
toionization calculations because those are aimed at pro- 
viding the emission lines without consideration of the 
shock heating in the first place, we envision a scenario 
where shock heating, despite the intense energy input 
from the supernova and stellar winds, is only a small 
fraction of the total energy output of the star formation. 
For our zeroth order model of star formation proceeding 
for 10 8 years, we find that t his condition is clearly sat- 
isfied (|Leitherer et al. I fl999). Thus we expect the emis- 
sion lines to be dominated by the ionizing energy output 
from massive stars and not by the shock heating. The 
clouds could therefore be thought of as the ones that have 
yet to be destroyed or are simply the surfaces of sheets 
and filaments that will get eventually or are now being 
run over by the blast waves into the ISM. 

The relationship for the a cceler ation of the 
cloud is given by IKlein et al. I (|1994 but see also 
INakamura et al. II2006I ). as^ 

m c ^ = -l/2C Dai « 2 A (6) 

where, m c is the mass of the cloud, ^ is the acceler- 
ation of the cloud, Cd ~ 1 is the drag coefficient, p^i is 
the post-shock density and is approximately, ~ 
(we are assuming strong shocks), the pre-shock density 
of the inter-cloud medium, and v c =| i^i — v c | is the rel- 
ative velocity of the shocked medium, where Vii is the 
mean velocity of the cloud, and A is the cross-sectional 
area. 

Does the star formation actually generate blast waves 
of sufficient energy to accelerate the clouds to the ob- 
served velocities? The total mass ejected by super- 
novae and stellar winds in the star formation process, 
Msf can be scaled as M = {3Msf, where is the 
fraction of the total ejected mass that is entrained or 
mass loaded in the winds. Simple conservation argu- 
ments suggest the terminal velocity of the outflow is 
Voo = (£:/M) - 5 « 2800(e//3)°- 5 km s" 1 , where e is the 
thermaliz ation efficiency of the starburst mechanical en- 
ergy (e.g., iMarcolini et al.ll2005l : IStrickland &: Heckman I 
2009). If we assume mass loading rates of 1-10, then the 
blast wave speed is about 1000-2000 km s _1 . Similarly, to 
first order the density of the blast wave is proportional to 
/g 3 / 2 which would suggest densities of-order 10~ 3 cm -3 
(e.g. . IMarcolini et al.ll2005t ). Such low densities and high 
velocities meet the general criteria of being a blast wave, 
which is highly supersonic for all of the densities in the 
ISM and should be very efficient in destroying clouds as 
it passes through the ISM. 

The modeling of the emission line gas in these galaxies 
suggest parameters for the clouds like density and column 
density, and assuming the ISM surrounding the pre-blast 
is like that typical of the ISM in local galaxies, implies 
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log ph=-3 to 0, log Pcioud=l to 3 cirT 3 , r cioud =0.1-10 
pc. Using these values for the clouds, the initial condi- 
tions for the ISM, and the blast wave formulation argued 
for previously as input to the simple cloud acceleration 
models, we find that the clouds could be accelerated up 
to 100 to a few 100 km s -1 before being destroyed. Thus 
it appears that at least in principle, the velocities of the 
outflows generated by intense star formation appear able 
to induce velocities like those observed. 

7.6. Bulk and Turbulent motions 

Ultimately, following the above reasoning, star forma- 
tion may give rise to the observed line widths through a 
combination of bulk and turbulent motion. This would 
also help resolve an apparent inconsistency in our above 
arguments. If turbulent motions dominate the observed 
velocity dispersions, we may expect the relationship be- 
tween the dispersion and star-formation intensity of the 
form, a—(eT,sFR) 1 ^ (ignoring the relationship between 
star-formation intensity and gas surface density which 
would tend to flatten this relationship). In Fig. [TJl 
we show this relationship based on Eqn. 3. To make 
this comparison with the data, we assume that the en- 
ergy of star-formation is dissipated entirely as turbu- 
lence, i.e., Ekin = Esf (see § 17.5. ip . assuming the 
same parameters used in § 17.11 and for two coupling ef- 
ficiencies between the energy injected by stars with the 
ISM. We have ignored the increase in the gas surface 
density with star-formation (the Schmidt-Kennicutt re- 
lation). While this relationship does explain the over- 
all trend in the data, data for individual galaxies are 
steeper than o~=(£Y,sfr) 1 ^ and in better agreement with 
a=(eY<sFR) 1 ^ 2 ■ It could be this mixture of bulk and tur- 
bulent velocities in the warm ionized gas steepens the 
relationship between the dispersion and star-formation 
intensity from a—(eT,sFR) 1 ^ 3 ■ But of course, as with the 
scaling a=(eT,sFR) 1 ^ 2 , this is too simplistic. One would 
need to consider the scale of energy injection, what is 
the true nature of the warm ionized gas, the gas den- 
sity, the dependence of the distribution of the gas phases 
with star- formation intensity, etc., in order to understand 
completely the underlying mechanisms exciting the gas. 

We emphasize that this analysis does not apply to all 
the emission line gas, only the highest surface bright- 
ness gas. In nearby starburst galaxies, very high velocity 
gas i s observed, up to several hundred to 1000 km s" 1 
(e.g. , iHeckman et~allll990l : iLehnert fc Heckman Ifl996ah . 
However, such gas is generally of low surface bright- 
ness, well b elow the detection limi t of the data presented 
here ( e .g., IHeckman et ah| |1990t ILehnert fe Heckman I 
Il996al ILehnert et al. I1999T) . As pointed out in 
e.g.. [Heckman et al.l ( 19901 ) and ILehnert fe Heckman I 
(1996a), the pressure in the emission line gas outside of 
the intense star-forming regions in local starburst galax- 
ies drops as roughly radius -2 and its surface brightness 
drops very rapidly as well. Only the nuclear regions with 
extremely high pressures reach surface brightnesses as 
observed in these distant galaxies and the molecular gas 
in such regions can shar e similar outflow velo cities as the 
warm ionized gas (e.g.. [Walter et al~l I2002L justifying 
our assumption that the kinematics of the warm ionized 
and molecular gas may be similar on the largest scales). 
We would therefore not expect to see the highest velocity 



gas in H a or in the high and low io nization lines observed 
(sec also iWang et al. l fl998L fl999T ) . 

Turbulence is thought to be driven by large scale bulk 
motions induced by energy injection - blast waves gen- 
erated by intense star-formation and other processes 
- over a wide range of scales, which then cascades 
into the denser gas and redistributes the energy over 
all scales in the med ium finally dissipating mainly on 
the smallest scales (Doung. Mac Low, fc Brvan I [2008; 
iMac Low fc Klessen Il2004 and references therein) . The 
nature of turbulence emphasizes the importance of both 
bulk and random motions and estimating energy dissi- 
pation rates based on turbulence arguments while esti- 
mating velocities based on bulk motion arguments are 
perhaps appropriate. Most likely, and in analogy with 
local starbursts, the warm ionized gas is probing the in- 
terface between outflowing gas and dense clouds in the 
ISM of these galaxies (especially at the highest surface 
brightnesses). It is therefore probing the interface where 
bulk motion and thermal energy is transfered to denser 
phases of gas through several mechanism such as ther- 
mal instabilities, drag against the background flows, col- 
lisions of cloud fragmen ts, gas cooling from warm HII 
to HI to denser H2 (see iGuillard et al. I [20091 . and ref- 
erences therein). Of course this requires efficient con- 
version of bulk moti ons into turbulen t energy which is 
apparently the case (|Mac Low Ill999h . The complexity 
of the processes that control the dynamics and distribu- 
tion of its various phases suggest that no simple scaling 
like o~=(eY>sFR) 1 ^ 2 can provide an ultimate understand- 
ing of the ISM. However, it does suggest that the ISM of 
these distant galaxies is controlled by the intense energy 
output of the star-formation within them. Not surpris- 
ing, but something that needed to be investigated and 
certainly needs further study. 

8. FURTHER IMPLICATIONS OF INTENSE 
STAR-FORMATION 

It appears that it is the intense star formation (or per- 
haps to some extent and in some cases the AGN) that is 
controlling the properties of the emission line gas within 
these galaxies. Given that we have found that the star 
formation has sufficient mechanical energy output in its 
own right to explain the characteristics of the emission 
line gas, it is unclear whether the gas is telling us much 
about the underlying mass dis tribution of these ga laxies 
or their origins. For example. iGenzel et a l. (2008) have 
argued that high central velocity dispersions require mass 
concentrations that are consistent with bulges. We can 
now explain this using the intense mechanical energy out- 
put of massive stars (with a non-negligable contribution 
from AGN in some cases like BzK-15504) and thus the 
velocity dispersions do not appear to be a unique tracer 
of the underlying gravitational potential. Similarly, it is 
plausible that variations in the physical conditions of the 
gas will lead to variation in the line ratios (see discus- 
sion in § 16.11 and Fig. [TT]) . Although these variations do 
not appear sufficient to significantly affect the metallic- 
ity estimates based on i ntegrated spectra (as argued by 
iBrinchmann et aT] [2008). they are sufficiently large that 
they may affect metallicity gradient estimates within in- 
dividual high redshift galaxies. 

Determining that the ISM of these distant galaxies ap- 
pears regulated by the mechanical output of the intense 
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star-formation, it is logical to investigate how this out- 
put might influence the nature of the star-formation it- 
self and our understanding of the dynamics of distant 
galaxies. We do not as yet have a comprehensive the- 
ory of star formation or a complete understanding of 
all the processes that may limit both the efficiency and 
the rate at which stars can form. Global star-formation 
laws or thresholds are those that describe the gross char- 
acteristics of star formation on large scales and over 
dynamical times of galaxies. Several possibilities for 
explaining the global characteristics of star for mation 
have been sugges ted: cloud-cloud collisions (e.g. JLarson I 
Il988t iTan I 20001) and the growth of gravitational per- 
turbations (jToomre II1964I). pre ssure and turbulenc e reg- 
ulated media (e.g.. isilk Ml997l 120011: lElmegreen 1 120021 : 
iMac Low fc Klessen 1120041) and others. Models such as 
cloud-cloud collisions and pressure and turbulence regu- 
lated ISM seem particularly appropriate in understand- 
ing the role of self-regulation in intense star-formation. 
However, of particular interest in distant galaxies, and 
something that we can directly comment on given our 
results, are the questions of the star-formation efficiency 
and how does intense star-formation affect the structural 
properties of distant galaxies. 

8.1. Efficiency of Star- Formation 

Recently, iDekel et al. I ([2009b) investigated whether 
or not cosmological gas accretion can sustain the star- 
formation rate of distant galaxies. Galaxies with high 
star formation rates as observed here require the star- 
formation rate to be close to the accretion rate of gas in 
this type of model. This, on the face of it, requires the gas 
infall time scale and the star-formation time scale to be 
close to equal which then implies that the star-formation 
efficiency is high. 

Our results suggest that the star-formation intensity, 
gas pressure, and velocity dispersions in z ~ 2 star- 
forming galaxies are similar or higher than those in 
the most intense starbursting regions at z = 0. The 
efficiency of star-formation in local spirals is about a 
few percent , but can re ach ~ 10% in starburst re- 
gions (e.g., lYoung et all 119961 : ISolomon fc Sage! 119881 : 
iSanders et all I1991D . From this, we would infer effi- 
ciencies of at least 10% in our z~2 galaxies. As men- 
tioned previously, we do not have a comprehensive the- 
ory of star formation. However, it has been suggested 
that there is a power-law dependence o f the star for- 
matio n rate on the total column density (jWang fc Silk I 
1994). Such a relationship implies that star- formation 
occurs on a dynamical ti me scale averaged over a large 
area. lElme green I (|2002f ) used this dependency to ar- 
gue that this would lead to a star-formation law of the 
form, SFI — esF^gasU, where csf is the star-formation 
efficiency and ui is a dynamical rate for the conversion 
of gas into dense , star- forming cloud cores. Interest- 
ingly, lElmegreen I ([2002D argue that this does not im- 
ply that gravitational forces are directly involved but 
that this dynamical rate of conversion is about equal 
to the turbulence crossing rate (which is also argued to 
be the inverse of the collapse time with modest over- 
pressures). Comparing the inverse of the relative cross- 
ing times of turbulence of our sample and local galaxies, 
(vturb/l)z~2/(vtur b/l)z~o ^2-5 for zwi,=10Q-2 00 km 
s _1 and 1~1 kpc (|Elmegreen fc E lmegreen 2006) in the 



distant galaxie s and ft MT .b=20— 25 km s _1 and 1~500 pc 
in local disks (|Epinat et al. I [2009] ). This suggests that 
the conversion efficiency of gas into dense star-forming 
cloud cores is higher in our sample of galaxies than for lo- 
cal disks. Thus the apparent efficiency of star-formation 
in our z^2 galaxies should be higher than in nearby 
disk galaxies, where by "apparent" we mean the prod- 
uct, esF^. But it is important to note that the star- 
formation efficiency itself, csf, can remain constant. On 
the other hand the efficiency can hardly exceed, say, 30%, 
since the star formation e fficiency at the scale of in divid- 
ual cores is about 50% ([Matzner fc McKee 1 12000ft and 
the efficienc y in molecular clo uds is necessarily some- 
what lower (|Elmegreen I l2002h . Our observations thus 
suggest "apparent" star-formation efficiencies between 10 
and 30%, not just in a nuclear region but over the whole 
disk, which is globally starbursting and highly turbulent. 

An increase in the apparent efficiency of star- formation 
can be understood with relation to the likely high sur- 
face density of gas in these galaxies. A high gas surface 
density is supported by both the intense star-formation 
(through the Schmidt-Kennicutt relation) and by the 
high pressures inferred from the optical emission line 
gas. If this efficiency is related to the inverse of the 
gas consumption time and the Schmidt-Kennicutt rela- 
tion, then the gas suface density divided by the star- 
formation intensity, the gas consumption time would be, 
Ytgas/YtsFR oc S~ Q ° S 5 - a decreasing function of the gas 
density. Thus finding a high apparent efficiency of star- 
formation is a natural consequence of having faster dy- 
namical processes such as faster collapse of gas due to 
high pressures, faster and greater turbulent compression, 
and stronger influence of self-gravity. Within this frame- 
work, it is easy to understand how the galaxies could 
be clumpy and highly unstable against rapid and intense 
star-formation. 

8.2. Clumpy Disks and Dynamical Mass Estimates 

Our results also h ave implications for the "cl umpy 

disk hypothesi s" ( Elmegreen fc Elmegreen I 120051 : 

lElmegreen et al.l 12007 ). whereby galaxy evolution is 
hypothesized to be driv en by internal clump formation 
(|Bournaud et al. I I2007D . Our results suggest that 
the clumpy disk model alone cannot explain the high 
dispersions observed and that an additional energy 
source is needed. Gravity can play a crucial role, 
generating the high masses o bserved in the clumps 
([Elmegreen fc Elmegreen 1 120051 ) . but we have argued 
that star-formation must be an important source of 
energy too. The intense star-formation is a mechanism 
whereby further clump formation m ay be stimulated 
by ma intaining the high dispersions. iBournaud et al.~l 
(2008) conclude that the observations of the "sky- 
walker" in the UDF (UDF6462) and that the peculiar 
velocities observed are consistent with a clumpy disk. 
Peculiar clump motions , typically around 50 km s" 1 
([Bournaud et al. 1 l2007l ). are likely masked by the 
turbulence and bulk velocities generated in the optical 
emission line gas by the intense star-formation. This 
implies that clump masses are difficult to infer from the 
dispersions of their optical emission line gas. Further- 
more, while the clumpy disk model is certainly viable 
for a subset of the observed galaxies, it is also likely that 
the mechanical output from the intense star-formation 
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within the clumps should not be ignored either and 
may have a profound influence on the evolution of the 
cl umps. 

iNoguchil (|1999f ) suggested that the bulges in local 
spiral galaxies form through the dynamical evolution 
of massive chimps of stars in high redshift galaxies. 
This m odel was further investigated bv lElmegreen et "all 
(2008a): large massive clumps migrate towards the disk 
center, and coalesce into a slowly-rotating bulge. The 
limited observations ava ilable may support this picture 
(|Elmegreen et al.~l l2009f ) . Could the high efficiency of 
star formation in our observations of z*~2 clumpy galaxies 
question this model? Vigorous energy injection from in- 
tense star formation may disrupt the clumps just like this 
process disrupts the star-forming clouds in low redshift 
galaxies. Actually, if the stars are forming in a higher 
density medium, which our observations may suggest, 
the clumps may be more tightly bound and more diffi- 
cult to destroy with the mechanical energy output from 
massive stars (even w ith strong feedback fro m intense 
star formation, e.g., iTasker fe Bryan 1 l2008f h There- 
fore, it is likely that the energy output affects mainly the 
lower- den sity, inter-clump g as, but not the most massive 
clumps. El megreen et al.l (|2008aD have independently 
argued that star-formation feedback could regulate the 
bulge growth, but does not prevent this bulge-forming 
mechanism from taking place in high redshift galaxies. 

Finally, we have argued that the emission line ratios 
for BzK-15504 are consistent with it harboring an AGN, 
and that this AGN is likely growing very rapidly. In this 
situation, it is difficult to know how much of the extended 
emission line gas is excited by star formation or the AGN 
and what drives the dynamics of the extended gas. BzK- 
15504, since it has an exquisite near-infrared data set, 
has been argued to be the archetypal gr owing disk at 
high redshift. In fact. lGenzel et all (|2008l ) have recently 
proposed that this galaxy, and some others, represent a 
phase of rapid gas accretion from their dark matter halo 
that is feeding its intense star formation. If the extended 
emission line gas is influenced by a powerful AGN, then 
this cannot be the case. The necessary gas masses and 
growth rates would be lower, perhaps by factors of few or 
more, than previous estimates if the AGN is increasing 
the overall surface brightness of the extended emission 
line gas. We do not think that the AGN is the only 
power source, but our data do not support the notion 
that this galaxy must have an extreme accretion rate or 
a rapid rotation suggestive of a massive disk. The ioniza- 
tion, surface brightness distribution, emission line distri- 
bution asymmetry, and line widths are all consistent with 
there being a kpc scale narrow line region surrounding a 
powerful QSO in this object. 

9. CONCLUSIONS 

Our analysis of integral field spectroscopic data with 
SINFONI on the ESO/VLT and HST/NICMOS imaging 
of our sample of z~2 galaxies have revealed a number 
of interesting properties and relationships. Our results 
suggest that the feedback cycle between the ISM and 
star formation (and perhaps, more speculatively, AGN) 
plays an important, perhaps crucial role in regulating 
star formation at high redshift. In support of this overall 
conclusion, we find: 

(1) These galaxies have surface brightnesses that are 



more than an order of magnitude greater than those of 
local disks. Our observations suggest that the ISM of 
these z~ 2 galaxies is starbursting over their whole area. 

(2) These galaxies have line ratios that suggest the pres- 
ence of an AGN, as in the case of BzK-15504, or high 
pressures and density bounded clouds. 

(3) Cosmological gas accretion models cannot provide 
sufficient energy if its energy dissipates as compressible 
turbulence. 

(4) A relationship between the star-formation intensity 
and the velocity dispersion of the emission line gas is 
found and it can be explained by a simple energy in- 
jection relation. At low velocity dispersions, self-gravity 
may play a role in generating turbulence. 

(5) The explanations for the relationships we observe im- 
ply that the Ha surface brightness and the distribution of 
velocity dispersion may tell us little about the underlying 
mass distribution. Thus bright Ha emission and broad 
lines are likely not evidence for mass concentrations such 
as bulges or rings within these galaxies. 

(6) Given the high pressures and velocity dispersions we 
observe, it is likely that the ISM (and perhaps star- 
formation) is turbulence an d pressure moderated (e.g., 

[19971 [20011 lElmegreerill200l [Mac Low fc KlesserTl 



Silk 


1997. 


2004 


)■ 



possi- 
2002; 



(7) Our results may imply apparent star-formation ef- 
ficiencies as high as 10% to 30% (for the stellar mass 
formed per gas mass unit per dynamical time). 

Through several lines of argument, we outlined a pic- 
ture in which self-gravity generates high turbulence in 
gaseous disks. Massive, large and gas-rich clumps then 
form, triggering star formation with high intensities 
(>0.1 M kpc -2 ) and apparent efficiencies (10-30%), 
with the "starbursting" interstellar medium spreading 
over scales of 10-20 kpc. At this stage of intense star 
formation, the turbulence and pressure moderated ISM 
likely regulates the st ar formation through several 
ble mechanisms (e.g. [Silk Ml997l l2001h lElmeg] 
iMac Low fc Klessen |[200l . The massive clumps how- 
ever likely survive and form bulge s following the mecha- 
nisms propose d bv INoguchil (Il999f) and discussed in fur- 
ther detail bv lElmegreen et al.1 (|2008al ). Such a picture 
has a myriad of implications for our understanding of the 
most rapid periods of star formation in distant galaxies 
and the ensemble properties of galaxies generally. 
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Intensely Star-Forming Galaxies at Redshift 2 



TABLE 1 
Properties of the High-z Galaxies 



Object 


Line 


z 


FWHM 


flux 


SB limit 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


O2343-BA610 


Ha 


2.2098±0.0018 


382± 8 


7.86±0.15 


4.2 


9.3 




[Nn] 


2.2104±0.0019 


4UU± 21 


3.23±0.15 








[Om] 


2.2090±0.0015 


3U1± 30 


0.67±0.06 








H/3 


2.2096±0.0019 


383± 02 


0.50±0.07 






ZC782941 


Ha 




305± 13 


5.23±0.21 


8.9 


7.2 




[Nil] 


2.1816±0.0017 


344± 75 


1.15±0. 22 






ZC782941'' 


Ha 


2.1814±0.0017 


346± 13 


3.05±0.11 


3.8 


4.2 




[Nil] 


2.1815±0.0018 


378± 42 


1.09±0.11 






Q2343-BX528 


Ha 


2.2684±0.0015 


293± 8 


2.28±0.05 


4.3 


7.5 




[Nil] 


2.2689±0.0017 


334± 63 


0.43±0.06 






BzK-15504 6 


Ha 


2.3816±0.0014 


260± 8 


6.02±0.16 


4.9 


5.8 




[Nil] 


2.3818±0.0017 


327± 38 


2.54±0.22 






BzK-15504 


Ha 


2.3819±0.0020 


400± 12 


7.58±0.20 


8.4 


9.3 




[Nn] 


2.3813±0.0032 


661± 95 


3.34±0.38 








[Om] 


2.3833±0.0024 


480± 15 


8.68±0.24 








H/3 


2.3841±0.0012 


182± 50 


0.73±0.17 






BzK-6397 


Ha 


1.5132±0.0011 


257± 8 


3.18±0.09 


2.3 


10.5 




[No] 


1.5135±0.0015 


406± 46 


1.08±0.11 






BzK-6004K 


Ha 


2.3865±0.0013 


245± 10 


4.28±0.16 


7.7 


8.7 




[Nil] 


2.3866±0.0015 


274± 22 


2.22±0.16 






Q2346-BX482 


Ha 


2.2563±0.0016 


312± 6 


4.18±0.07 


4.4 


9.0 




[No] 


2.2564±0.0016 


313± 40 


0.64±0.07 








[Om] 


2.2569±0.0012 


212± 11 


1.62±0.07 








H/3 


2.2589±0.0011 


166 


0.22 






K20-ID7 


Ha 


2.2234±0.0016 


313± 9 


4.02±0.10 


5.3 


9.3 




[Nn] 


2.2242±0.0014 


282± 50 


0.79±0.11 






K20-ID8 


Ha 


2.2231±0.0014 


267± 10 


2.51±0.09 


7.2 


7.5 




[Nn] 


2.2238±0.0014 


264± 50 


0.75±0.10 






SSA22a-MD41 


Ha 


2.1632±0.0015 


298± 12 


3.71±0.13 


5.2 


8.3 




[Nn] 


2.1627±0.0018 


382± 89 


0.72±0.16 








[Om] 


2.1704±0.0016 


323± 14 


2.29±0.09 








H/3 


2.1704±0.0009 


96± 19 


0.44±0.07 






Q2343-BX389 


Ha 


2.1716±0.0024 


525± 18 


5.44±0.17 


8.2 


8.5 




[Nn] 


2.1732±0.0031 


681±112 


1.65±0.22 








[Om] 


2.1713±0.0019 


393± 24 


1.65±0.09 








H/3 


2.1708±0.0010 


141± 28 


0.38±0.06 







Note. — Column (1) - Object designation. a implies 250 mas pixel scale plus 
adoptive optics, b implies 100 mas pixel scale plus adaptive optics. No indication implies 
that the data were taken without the benefit of AO and at 250 mas pixel^ 1 . Column (2) 

- Line identification. Column (3) - Redshift of the line in the integrated spectrum. By 
integrated spectrum, we mean that the sum of the flux from each pixel with a signal-to- 
noisc greater than or equal to 3 in each data cube for Ha (see column 6 for this limiting 
value for Ha). All of the sums for the other emission lines are over the same aperture 
as for Ha. Column (4) - Full width at half maximum of the integrated spectrum of 
each galaxy, corrected for instrumental resolution and is in units of km s _1 . Column (5) 

— Line flux of the integrated spectrum in units of 10 — 16 erg s _1 cm -2 . Column (6) - 
Surface brightness detection limit in units of 10 -19 erg s _1 cm -2 and defined at a signal 
to noise ratio, S/N?s3. This is for a pixel that has been averaged over 3 pixels x 3 pixels. 
Column (7) - Isophotal radius defined as Area — it r 2 where the Area is defined as the 

2 

projected area on the sky above a significance of 3 in the data and is in units of kpc. 
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TABLE 2 
H/3 Line Fluxes and Extinctions 



Source 
(1) 


ft 


fffa/fff/3 

(3) 


A H/3 
(4) 


corr Ha 
(5) 


Q2343-BX610 


0.61±0.08 


8.1 


2.8 


4.6 


Q2343-BX389 


0.37±0.05 


4.1 


1.0 


1.7 


SSA22a-MD41 


0.35±0.05 


7.2 


2.5 


3.9 


BzK-15504 


0.52±0.08 


7.1 


2.4 


3.8 



Note. — Column (1) - Source designation. Column (2) - 
Measured line flux of H/3 in units of 10 -16 erg s _1 cm -2 for 
the integrated spectrum. In this case, we mean that the spec- 
trum is integrated over the region where the S/N is greater 
than 3 in H/3. Column (3) - Line ratio {hoc/{h& f° r the inte- 
grated spectrum as defined for column 2 in this table. Column 
(4) — The extinction in the H/3 using the Galactic extinction 
law. Column (5) — The multiplicative factor required to cor- 
rect Ha for extinction assuming the Galactic extinction law. 



TABLE 3 

[SII] Emission Line Fluxes and Electron Densities 



Source 


fA6716 


fA6731 


fA6716/fA6731 


n e 




(1) 


(2) 


(3) 


(4) 


(5) 




Q2343-BX389 


5.9±0.1 


7.5±0.1 


0.8±0.1 


1200 


+700 

-400 


Q2343-BX610 


1.7±0.1 


1.5±0.1 


l.liO.l 


400" 


-700 
-300 


Q2347-BX482 


2.4±0.1 


3.1±0.1 


0.8±0.1 


1200 


+ 700 
-400 


BzK-15504 


2.9±0.1 


2.8±0.1 


1.1±0.1 


4ooi™:: 


BzK-6397 


l.liO.l 


0.9±0.1 


1.2±0.1 


260: 


-150 
-+20 



Note. — Column (1) - Source designation. Column (2) - Mea- 
sured line flux of [Sll] A6716 in erg s~ 1 cm -2 . Column (3) - Measured 
line flux of [Sn]A6731 in erg s _1 cm -2 . Column (4) - Line ratio 
F(6716)/F(6731). Column (5) - Electron density corresponding to 
R±lcr for T= 10 4 K in cm" 3 . 



